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ABSTRACT 

The increasing richness of exoplanet observations has motivated a variety of three-dimensional atmospheric 
circulation models of these planets. Under strongly irradiated conditions, models of tidally locked, short-period 
planets (both hot Jupiters and terrestrial planets) tend to exhibit a circulation dominated by a fast eastward, or 
"superrotating," jet stream at the equator. When the radiative and advection time scales are comparable, this 
phenomenon can cause the hottest regions to be displaced eastward from the substellar point by tens of de- 
grees longitude. Such an offset has been subsequently observed on HD 189733b, supporting the possibility 
of equatorial jets on short-period exoplanets. Despite its relevance, however, the dynamical mechanisms re- 
sponsible for generating the equatorial superrotation in such models have not been identified. Here, we show 
that the equatorial jet results from the interaction of the mean flow with standing Rossby waves induced by 
the day-night thermal forcing. The strong longitudinal variations in radiative heating — namely intense day- 
side heating and nightside cooling — trigger the formation of standing, planetary- scale equatorial Rossby and 
Kelvin waves. The Rossby waves develop phase tilts that pump eastward momentum from high latitudes to the 
equator, thereby inducing equatorial superrotation. We present an analytic theory demonstrating this mecha- 
nism and explore its properties in a hierarchy of one-layer (shallow- water) calculations and fully 3D models. 
The wave-mean-flow interaction produces an equatorial jet whose latitudinal width is comparable to that of the 
Rossby waves, namely the equatorial Rossby deformation radius modified by radiative and frictional effects. 
For conditions typical of synchronously rotating hot Jupiters, this length is comparable to a planetary radius, 
explaining the broad scale of the equatorial jet obtained in most hot Jupiter models. Our theory illuminates the 
dependence of the equatorial jet speed on forcing amplitude, strength of friction, and other parameters, as well 
as the conditions under which jets can form at all. 

Subject headings: hydrodynamics - methods: analytical - methods: numerical - planets and satellites: atmo- 
spheres - planets and satellites: individual (HD 189733b) - waves 



1. INTRODUCTION 

The past few years have witnessed major strides in our ef- 
forts to understand the atmospheric circulation of short-period 
exoplanets — both gas giants (hot Jupiters) and smaller terres- 
trial planets. Infrared photometry, spectra, and light curves 
from the Spitzer and Hubble Space Telescopes now provide 
constraints on the three-dimensional temperature structure of 
several hot Jupiters, which hint at a vigorous atmospheric 
circulatio n on these bodi es (e.g., Knutson et al. 2007, 2009| 
Charbonneau et al.|2008[ [Harrington et al. ^2006^ Cowan et alT] 
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servations have motivated a growing effort to model the at- 
mospheric circulation on these objects: to date, many three- 
dimensional atmosph eric circulation models of hot Jupiters 
have bee n pubHshe d ([Showman & Guil lotll2002( [ Cooper & ] 
Showman 2005, 2006; Showman et al. 2008, 2009; Dobbs- 
Dixon & Lin, 2008; Menou & Rauscher 2009; Rauscher & 
Menou||2010T|D obbs-Dixon et al.||2010t [T hrastarson & Cho 
2010t [Lewis 'et al.|2010nPem a et al.'2010; He ng et al.|2010^ . 
These models have emiphasized synchronously rotating hot 
Jupiters in circular, approximately 2-5 -day orbits. 
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Just as the last decade witnessed the first characteriza- 
tion of hot Jupiters, the next decade will see a shift toward 
characterizing "super Earths" (planets of 1-10 Earth masses) 
and terrestrial planets. To date, roughly 30 super Earths 
have been discovered, inc luding several that transit their host 



stars jCharbonneau et al . 2009; Leger et al. 2009; Batalha| 
et al.||2011| ) with hundreds of additional cand idates recently 



annou nced from the NASA Kepler mission (Borucki et al. 
2011| ). Attempts to observationally characterize their atmo 



spheres have already begun (Bean et al. 2010). In anticipa- 
tion of this vanguard, several three-dimensional circulation 
models of tidally loc ked, short-period terrestrial exoplanets 
have been published (|Joshi et al.|1997l|Joshi|2003||Merlis & 



Schneider|201QllHeng & Vogt|2010l) . 

Intriguingly, tne flows in most of these three-dimensional 
models — both hot Jupiters and terrestrial planets — develop 
a fast eastward, or superrotating, jet stream at the equator, 
with westward flow typically occurring at deeper levels and/or 
higher latitudes. In hot- Jupiter models, the superrotating jet 
extends from the equator to latitudes of typically 20-60° and 
is perhaps the dominant dynamical feature of the modeled 
flows. In some cases (depending on the strength of the im- 
posed stellar heating and other factors), this jet causes an 
eastward displacement of the hottest r egions from the sub- 
stellar point by '^10° to 60° longitude. [Showman & Guillot 



( 2002| ) first predicted this feature and suggested that, if it ex 



isted on hot Jupiters, it would have important implications for 
infrared spectra and light curves. This prediction has been 
confirmed in Spitzer infrar ed observations of HD 189733b 
( [Knutson et al.||2007[ [2009| ), suggesting that this planet may 
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indeed exhibit a superrotating jet. 

However, despite the ubiquity of equatorial superrota- 
tion in three-dimensional models of synchronously rotating 
exoplanets — and its relevance for observations — the mecha- 
nisms that produce this superrotation have yet to be identified. 
As demonstrated in a theorem due to Hidej(^1969 ), such super- 
rotation cannot result from atmospheric circulations that are 
longitudinally symmetric or that conserve angular momentum 
per mass about the planetary rotation axis. The equatorial at- 
mosphere is the region of the planet farthest from the plane- 
tary rotation axis, and a superrotating equatorial jet therefore 
corresponds to a local maximum in the angular momentum 
per mass about the planetary rotation axis. Thus, any angular- 
momentum-conserving circulation that moves air to the equa- 
torial atmosphere from higher latitudes or deeper levels tends 
to produce westward equatorial flow. Equivalently, Coriolis 
forces always induce westward acceleration for air moving 
equatorward or upward, so an eastward equatorial jet can- 
not result from Coriolis forces acting on air that moves into 
the equatorial atmosphere from higher latitudes or deeper lev- 
els. In Earth's equatorial troposphere, for example, the flow 
is westward, which results from the tropospheric Hadley cell 
on Earth (a regime where a mean overturning circulatio n and 
its Co riolis accelerations plays the defining role (Held & Hou| 
1980| )). To maintain equatorial superrotation, a mechanism is 
needed that pumps angular momentum per mass from regions 
where it is small (outside the jet) to regions where it is large 
(within the jet) — a so-called "up-gradient" momentum trans- 
port. According to Hide's theorem, this transport can only be 
accomplished by waves or eddies. 

Equatorial superrotation exists in several atmospheres of 
the Solar System — the equatorial atmospheres of Venus, Ti- 
tan, Jupiter, and Saturn all superrotate. Even localized lay- 
ers within Earth's equatorial stratosphere exhibit superrota- 
tion, part of a phenome non called the '' Quasi-Biennial Os- 
cillation" or QBO ( An drews et al.lp^'STl ). The mechanisms 
for driving the equatorial superrotation on these planets are 
diverse. Possible mechanisms include eddy transport associ- 
ated with baroclinic instabilities, barotropic instabilities, tur- 
bulence, and the absorption/radiation of various types of at- 
mospheric waves (e.g., Williams 2003b a; |Lian & Showman 
2008, 201()l|Schneider & Liu 2009, Del Genio et al.|1993 ;.DeI 



tion reads 



Genio & Zhou| 19961 [Andrews et al.|1987t|Mitchen & Valhs 

loTol T 

Here, we demonstrate how the equatorial superrotation in 
three-dimensional models of synchronously rotating, short- 
period exoplanets can result from the existence of stand- 
ing, planetary- scale Rossby waves; such waves are naturally 
excited by the longitudinally dependent heating patterns — 
dayside heating and nightside cooling — that accompany the 
photospheric regions of short-period, synchronously rotating 
planets. Section 2 provides background. In Section 3, we 
present an analytic theory demonstrating the mechanism, and 
we systematically explore its behavior in idealized, nonlin- 
ear one-layer models. In Section 4 we extend the analysis to 
three-dimensional circulation models. In Section 5, we sum- 
marize and discuss implications. 

2. BACKGROUND THEORY 

The ability of Rossby waves to accelerate jets can be 
schematically illustrated using the two-dimensional non- 
divergent barotropic vorticity equation, which is the simplest 
model for the global-scale flow of a pla netary atmosphere (see 
discussion in Vallis^2006{ [Showman et al.^2010) . The equa- 



^(C+/) 

dt 



= F 



(1) 



where ^ = k • V x v is the relative vorticity, v is the hori- 
zontal wind velocity, k is the local upward unit vector, / = 
21] sin is the Coriolis parameter, ^} is the planetary rota- 
tion rate (In over the rotation period), (j) is latitude, and 
d / dt = d / dt -\-Y ■ \/ is the material derivative (i.e., the deriva- 
tive following the flow). The equation states that individual 
fluid parcels conserve the absolute vorticity, C + /' save for 
vorticity sources/sinks, which are represented by the term F. 
The equation can equivalently be written 



dC 
dt 



(2) 



where v is the meridional (northward) wind speed and f3 = 
df/dy is the gradient of the Coriolis parameter with north- 
ward distance y. Because the flow in this simple model is 
horizontally non-divergent, we can define a streamfunction 
such that u = -d^l) / dy and v = dil)/ dx, where x is the eastward 
coordinate and u is the zonal (eastward) wind speed. This al- 
lows the equation to be written as a function of one variable. 

For purposes of illustration, consider the linearized version 
of Eq. (|2]) with no sources and sinks. The solutions to this 
linearized, unforced equation are Rossby waves. For sim- 
plicity, we consider Cartesian geometry with constant /3, rep- 
resenting a local region on the sphere. Decomposing vari- 
ables into zonal means (denoted with overbars) and deviations 
therefrom (denoted with primes), and assuming that the mean 
flow is zero, leads to a solution given by ip^ = iljexp[i(kx-\-ly)], 
where / is the imaginary number and k and / are the zonal and 
meridional wavenumbers. The dispersion relation is 



(3k 



(3) 



These waves propagate meridionally with a group velocity 
given by 

dij _ ipkl 
~dl~W+t^' 



A si mple argument, first c learly pre sented by |Thompson| 
( [T97T] ) and reviewed in Heidi ( |20Q0l ) and |VaIIisl (|2006'), shows 
how these waves can produce an east- west acceleration of the 
zonal-mean flow. The latitudinal transport of eastward eddy 
momentum per mass is u'v' , where u' and v' are the deviations 
of the zonal and meridional winds from their zonal means, 
respectively, and the overbar denotes a zonal average. Using 
the solutions for the wave-induced zonal and meridional wind, 
u' = -ilil; exp[i(kx + ly)] and V = ikiljexp[i(kx + ly)], yields a 
momentum flux 



u'v' = --4)^kl. 



(5) 



Since the group velocity must point away from the region of 
wave generation (which we call the "wave source"), and since 
(3 is positive, we must have kl > north of the source and 
kl < south of the source. Thus, north of the source, u'v' 
is negative, implying southward transport of eastward mo- 
mentum. But south of the source, u'v' is positive, implying 
northward transport of eastward momentum. Rossby waves 
therefore transport eastward momentum into the wave source 
region. An eastward acceleration must therefore occur in the 



Equatorial superrotation on hot exoplanets 



3 



wave source region and a westward acceleration must occur 
in the region of wave breaking or dissipation. This can lead to 
the formation of zonal (east- west) jet streams r] 

Rossby waves correspond to latitudinal oscillations in sur- 
faces of constant potential vorticit}]^ thus, any process that 
triggers such oscillations at large scales will tend to excite 
Rossby waves. In Earth's atmosphere, one of the predomi- 
nant sources is baroclinic instability, which occurs in the mid- 
latitude troposphere where latitudinal temperature gradients 
are large. Spatially varying tropospheric heating and cooling 
(e.g., due to land- sea contrasts) or flow over topography also 
perturb the potential vorticity contours and can therefore trig- 
ger Rossby waves. In the atmospheres of tidally locked, hot 
exoplanets, on the other hand, the day-night heating pattern 
constitutes the overriding dynamical forcing. For such plan- 
ets, we expect this heating/cooling pattern to trigger Rossby 
waves at low latitudes (Fig.[T]). 

The above theory is for free waves. Consider now the ex- 
tension to an atmosphere forced by vorticity sources/sinks and 
damped by frictional drag. The zonal-mean zonal momentum 
equation of the barotropic system reads 



of the so-called "pseudomomentum" (I Vallis|2Q06 p. 493): 



du 
dt 



dy 



u 



(6) 



where overbars denote zonal means and primes denote devia- 
tions therefrom. The equation states that accelerations of the 
zonal-mean zonal flow result from convergences of the latitu- 
dinal eddy momentum flux and from drag, which we have 
parameterized as a term that relaxes the zonal-mean zonal 
wind toward zero over a drag time constant Tdrag. The rela- 
tionship between the eddy acceleration in Eq. ([6]) and the vor- 
ticity sources/sinks can be made in two steps. First, we note 
that the definition of vorticity implies that vX' = -d{u'v')/dy. 
Second, we multiply the linearized version of Eq. ^ by C' 
and zonally average. This leads to an equation for the budget 



^ The dynamical picture outlined above is not limited to small-amplitude 
di sturb ances, as can be shown with a simple argument described for example 
in "Held ( 2000) and Vallis ( 2006 1. Imagine an initially undisturbed latitude, 
where the absolute-vorticity contour initially aligns with the latitude circle, 
and suppose a disturbance — of any amplitude — propagates into that latitude 
from elsewhere. The disturbance will perturb the absolute vorticity contours, 
causing northward transport of air in some regions and southward transport in 
others. Because absolute vorticity generally increases northward, the north- 
ward advection carries with it air of low absolute vorticity, whereas the south- 
ward advection carries with it air of high absolute vorticity. Thus, this process 
will generally cause a southward flux of absolute vorticity, thereby decreas- 
ing the areal integral of the absolute vorticity over the polar cap bounded 
by the latitude circle in question. By Stokes' theorem, this implies that the 
zonal-mean zonal wind decelerates (i.e., accelerates westward) because of 
this vorticity flux. In absence of dissipative processes, this deceleration would 
reverse if the disturbance exited the region. However, when mixing occurs 
(e.g., if the wave breaks), or if the disturbance is damped before air parcels 
can return to their original latitudes, then the areal integral of the vorticity 
inside the latitude circle has been irreversibly decreased, and the westward 
impulse cannot be undone. Thus, we again recover the result that westward 
acceleration occurs in the region of wave dissipation; if momentum is con- 
served, eastward acceleration would then occur in the wave-source region. 

^ Potential vorticity is a quantity related to vorticity that is conserved in 
adiabatic, frictionless, stratified flow. For the barotropic system it is simply 
the absolute vorticity C + /; for the shallow-water system it is absolute vor- 
ticity over layer thickness ((^ + f)/h, and for a three-dimensional stratified 
atmosphere is it given by p~^(V x \ + 2fl) • V6, where p is density, fl is the 
planetary rotation vector, and is potential temperature. For discu ssion of 
the conservation of potential vorticity and its uses in dynamics, see |Pedlosky| 
fl987 > or VaUis ( 2006) . 



dA 
dt 



C'F' 



2(/5-p) 



(7) 



For the two-dimensional non-divergent model, A = - 
d^ujdy^y^C!^ 12 is the pseudomomentum, which is a mea- 
sure of wave activity. By combining Eqs. ([6])-([7]) and suppos- 
ing that the wave amplitudes and zonal-mean zonal wind are 
statistically steady, i.e., dA/dt ^ and du/dt ^ 0, we obtain 



(8) 



This equation relates the vorticity sources/sinks and drag to 
the zonal-mean zonal wind, u. When eddy sources/sinks of 
relative vorticity on average exhibit the same sign as the vor- 
ticity itself (i.e. C,'F' > 0), the eddy acceleration is eastward, 
and in steady state results in an eastward zonal-mean zonal 
wind. When sources/sink of relative vorticity tend to exhibit 
the opposite sign as the vorticity (C'F' < 0), the eddy accel- 
eration is westward, and in steady state results in westward 
zonal-mean zonal wind|^ In analogy with the free solutions, 
this behavior is typically interpreted in terms of the genera- 
tion, latitudinal propagation, and dissipation of Rossby waves. 

This mechanism is thought to be responsible for the 
eddy-driven jet streams (and the associated eastward surface 
winds) in Earth's midlatitudes: baroclinic instability gener- 
ates Rossby waves that radiate away from the midlatitudes, 
causing eastward eddy acceleration there and leading to east - 
ward surface flow in steady state ( |Held||2Q00t [Vallis ||20Q6"1 ). 
At the equator. Earth's troposphere does not superrotate; nev- 
ertheless, idealized Earth general circulation models (GCMs) 
have shown that the presence of strong zonally varying heat- 
ing and cooling in the tropics can cause equatori al superr ota- 
tion to emerge ([S uarez & Duffy 1992; S aravanan 1993 ; Krau-| 
cunas & Hartmann|2005 ; Norton 2006 ). In analogy with tEe 
theory described above, iHeld (1999) suggested heuristically 
that the superrotation in these models results from the gener- 
ation and poleward propagation of Rossby waves by the trop- 
ical heating sources. 

Still, application of this barotropic theory to tidally locked 
exoplanets is problematic. Most models of atmospheric 
circulation on synchronously rotating, zero-eccentricity hot 
Jupiters exhibit relatively steady circulation patterns whose 
velocity and temperature patterns are approximately synmiet- 



ric about t he equator ( Showman & GuiHot| |2002^> 'Coop er & 
Showman|r2 005, 20 061 [Showman et al.||200 8, 2009; DoW£ 
Dixon & Lin||2008nRauscher & Menou| |2010). In such 



models, the relative vorticity is approximately antisymmetric 
about, and zero at, the equator. Eq. ^ predicts that u = at 
the equator for this situation. Thus, while attractive, this the- 
ory fails to explain the equatorial superrotation in these hot 
Jupiter models. 

There are other difficulties. First, the theory presented here 
assumes the flow is barotropic — i.e. that it can be described 
by Eq. ([T]) — and thereby ignores the role of finite Rossby 
deformation radius in shaping the wave properties. Second, 
vertical motions normally accompany regions of heating and 
cooling, and these vertical motions lead to nonzero horizontal 
divergence. The low-latitude flow is thus inherently divergent 

^ These arguments assume that /S-d^u/dy^ > 0, which is generally the 
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within an equator ial waveguide and cannot propagate away 




Fig. 1 . — Illustration of the dynamical mechanism for generating equatorial 
superrotation on tidally locked short-period exoplanets, including hot Jupiters 
and super Earths. The intense day-night heating gradient generates standing, 
planetary- scale Rossby and Kelvin waves. These waves develop a structure 
with velocities tilting northwest-to-southeast in the northern hemisphere and 
southwest-to-northeast in the southern hemisphere (yellow and red ovals). 
In turn, these patterns transport eddy momentum from high latitudes to the 
equator (dashed arrows). Equatorial superrotation therefore emerges (thick, 
right-pointing arrow). 



in the presence of heating/cooHng, in conflict with the stated 
assumptions. Moreover, the theory cannot account in any rig- 
orous way for the generation of the Rossby waves by thermal 
forcing. The above equations include no thermodynamics and 
only the effect of such forcing in producing relative vorticity 
is represented. Away from the equator, one expects that heat- 
ing produces rising motion and horizontal divergence aloft, 
which generates anticylonic relative vorticity by the action of 
Coriolis forces. At the equator, however, this vorticity source 
is small, leaving unclear the applicability of this barotropic 
theory in producing eastward flow at the equatorij When 
the geostrophic assumption is relaxed and finite deformation 
radius is included, analytic solutions of freely propagating 
equatorial waves show that such waves tend to be trapped 

^ Within the context of barotropic theory, one can relax the nondivergent 
assumption by resolving the horizontal velocity into rotational and diver- 
gent cornponents, specifying the horizontal divergence field, and then solv- 
ing Eq. I?) for the rotational component of the flow (see Sardeshmukh & 
|Hoskins TP 85 1988 1. In this context, the specified divergence field repre- 
sents the spatial field of heating and cooling. The vorticity source in Eq. j2| 
would then be F = -((" -I- /)V • v^^ - v^^ • V(C + /), where v^^ is the specified 
divergent component of the flow field and the velocity and vorticity on the 
lefthand side of Ecl misrepresent only the rotational component. One can 
then rework Eqs. jS}-^ under the usual barotropic assumptions that verti- 
cal momentum/vorticity transport and vortex tilting are negligible. However, 
such a divergent barotropic model suffers from the same failing as the nondi- 
vergent version: when the pattern of heating/cooling and flowfield are mirror 
symmetric about the equator — as in most hot- Jupiter models — the vorticity 
source F is zero at the equator, leading to a Rossby- wave source d^'F' that 
is likewise zero at the equator. In the presence of frictional drag, one again 
obtains the result that w = at the equator. Thus, even such an extended 
barotropic treatment is insufficient to explain the equatorial superrotation in 
hot- Jupiter models. 



from the equator ( |Holton||2004t [Andrews et"ar||1987| ); such 
wave solutions involve no net nieridional (north- south) mo- 
mentum transports, thus raising the question of whether the 
above mechanism is viable at the equator in the presence of 
finite deformation radius. Finally, for hot Jupiters, the Rossby 
waves are expected to be global in scale, and it is not clear a 
priori whether there is room for them to propagate poleward 
from the equatorial regions. 

We present a sequence of models in the following sections 
that overcome these obstacles and provide a theoretical foun- 
dation for understanding equatorial superrotation on tidally 
locked exoplanets. 

3. SHALLOW- WATER MODEL OF EQUATORIAL 
SUPERROTATION 

Full GCM solutions, although useful, involve so many in- 
teracting processes that it is often difficult to cleanly identify 
specific dynamical m echanisms from such solutions (see, e.g., 
Showman et al.| !2Q10). Simplified models therefore play an 
important role in atmospheric dynamics. Here, we present a 
highly idealized model intended to capture the mechanism in 
the simplest possible context. 

We adopt a two-layer model, with constant densities in each 
layer, where the upper layer represents the meteorologically 
active atmosphere and the lower layer represents the quiescent 
deep atmosphere and interior. In the limit where the lower 
layer becomes infinitely deep and the lower-layer winds and 
pressure gradients remain steady in time (which requires the 
upper layer to be isostatically balanced), this two-layer sys- 
tem reduced to th e shallow-wa ter equations for the flow in the 
upper layer (e.g., IValhs 12006) p. 129-130): 



— +gVh+fkx^ = R 

at '?~drag 



dh /Zeq(A,(/))-/z 

-- + V • {yh) = 

at Trad 



Q 



(9) 



(10) 



where v(A,(/),0 is horizontal velocity, /z(A, (/),/_) is the upper 
layer thickness, t time, g is the (reduced) gravit}[j f = 2Q sin (j) 
is Coriolis parameter, k is the upward unit vector, Q is plane- 
tary rotation rate, and (j) is latitude. Again, d/dt = d/dt+\ -V 
is the material derivative. The boundary between the layers 
represents an atmospheric isentrope, across which mass is ex- 
changed in the presence of heating or cooling. Heating and 
cooling are therefore represented in the shallow-water system 
using mass sources and sinks, represented here as a New- 
tonian relaxation of the height toward a specified radiative- 
equilibrium height /Zeq — thick on the dayside and thin on the 
nightside — over a radiative time scale Trad- The momentum 
equations ^ include drag with a timescale Tdrag, which could 
represent trie potenti ally impo rtant effects of magnetohydro- 
dyn amic friction (Perna et al.'2010 ), vertical turbulent mixing 
(e.g. |Li & Goodm an 2010 ), or momentum t ransport by b reak- 
ing gravity waves (Lindz en,1981y,Watkins & Cho|2010| ). 

The term R in Eq. ([9]) represents the effect on the upper 
layer of momentum advection from the lower layer, and takes 

^ g in Eq. jo} is the actual gravity times the fractional density difference 
between the layers, and is therefore called the "reduced gravity." When in- 
terpreting this shallow-water model in the context of a three-dimensional at- 
mosphere, this density difference should be interpreted as (for example) the 
fractional change in potential temperature across a scale height. This is of 
order unity for a hot Jupiter. 
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the same form as in [Shell & Heidi ( |2QQ4| ) and [Showman "&] 
[Polvani[ ( [2QT0l ): 



0, 



(11) 



where A and cj) are longitude and latitude. Air moving out 
of the upper layer (2 < 0) does not locally affect the up- 
per layer's specific angular momentum or wind speed, hence 
R = for that case. But air transported into the upper layer 
carries lower-layer momentum with it and thus alters the local 
specific angular momentum and zonal wind in the upper layer. 
For the simplest case where the lower layer winds are assumed 
to be zero, this process preserves the column-integrated yh of 
the upper layer, leading to the expression in ([TT]) for 2 > 0. 
Importantly, the expression for R follows directly from the 
momentum budget and contains no free parameters n 

The relative roles of the two terms on the right side of the 
momentum equation can be clarified by rewriting it as fol- 
lows: 



dy 
dt 



+gV/z+/k X V = -V 



'^rad 



-h 



(12) 

where H(/zeq - h) is the Heaviside step function, defined as 1 
when hQq-h> and otherwise. Dynamically, it should be 
clear that both terms on the right side play a role analogous to 
drag; one can define the entire quantity in square braces as one 
over an effective drag time constant. Still, the second term (R) 
is spatially heterogeneous and only exists in regions of heat- 
ing, and we will show that its effect on the zonal-mean flow 
is qualitatively different than that of the first term (frictional 
drag). For a strongly irradiated hot Jupiter, we might expect 
(/Zeq -/?)//? ~ 0.01-1, and if so then the first term would dom- 
inate if Tdrag <C Trad whcrcas the second term would dominate 

if '^drag '^rad • 

We present linear, analytic solutions and fully nonlinear, 
numerically determined solutions of Eqs. ([9])-(pT|) in the next 
two subsections. 

3.1. Linear solutions 

To enable analytic solutions, we solve Eqs. ([9|)-(pT]) in 
Cartesian geometry assuming that the Coriolis parameter can 
be approximated as / = Py, where y is northward distance 
from the equator and (3 (the gradient of Coriolis parameter 
with northward distance) is assumed constant. This approxi- 
mation, called the "equatorial /3-plan e," i s strictly valid only 
at low latitudes, but we will see in |3.2[ that the qualitative 
features of these solutions are recovered by the full solutions 
in spherical geometry. 

We now linearize Eqs. ([9|)-(pT]) about a state of rest. By def- 
inition, all the terms in the linearized equations have magni- 
tudes that scale with the forcing amplitude. Note that the term 
R involves the product of the velocity and forcing amplitude, 
and therefore is quadratic in the forcing amplitude and does 
not appear in the linearized equations. (We will come back 

^ The condition that air moving out of the upper layer induces no change in 
the upper layer's specific momentum requires that such momentum advection 
cause no accelerations, implying that R = for Q < 0. To derive the expres- 
sion for the case Q > 0, add v times the continuity equation to h times the 
momentum equation, thus yielding an equation for the time rate of change of 
\h. Terms involving heating/cooling are hR + \Q. In the special case where 
the lower-layer winds are zero, mass transport into the upper layer does not 
change the column-integrated horizontal wind, \h, of the upper layer. This 
implies that hR + \Q = 0, thus yielding Eq. |7T|. 



to it when evaluating the implications of the solutions for the 
zonal-mean zonal wind). The linearized equations read 



du dri ^ u 
^+^^~^->^^" 

Ot ox Td 

dv dri ^ 



drj 

'dt 



+H 



du dv 
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'^drag 
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'^rad 



(13) 
(14) 

(15) 



where x is eastward distance and r] is the deviation of the 
thickness from its constant reference value //, such that h = 
H + T]. The quantity S = {h^c^- H) / r^^^ is the forcing, which 
can also be expressed as 5" = T^eq/'H-ad, where r/eq = hQq-H is 
the deviation of the radiative-equilibrium height from H. 

The free solutions to these equations (i.e., when the right- 
hand sides are set to zero) are the well-known e quatori- 
ally tr apped wave modes, described for example in [Holton 



( [2604[ ) and Andrews et al (1987 ). Given the intense heat- 
ing and cooling experienced by hot, tidally locked exoplan- 
ets, however, we seek solutions to the forced problem. Most 
three-dimensional dynamical models of hot Jupiters exhibit 



relatively steady circulation patterns (Showman & Guillot 
2002; Cooper & Showman 2005, 2006; Dobbs-Dixon & Lin 
2008; Dobbs-Dixon et al. 2010; Show man et al. 2008, 2009f 
Rauscher & Me nou 2010), and so we seek steady solutions in 
the presence of forcing and damping. 

We nondimensionalize Eqs. (p3])-(p3]) with a length scale 
L = {y/gH //3)^/^, a velocity scale U = ^/gH, and a time scale 
T = (y/gUpy^^^, which correspond respectively to the equa- 
torial Rossby deformation radius, the gravity wave speed, and 
the time for a gravity wave to cross a deformation radius in 
the shallow-water system. The thickness is nondimensional- 
ized with H, the drag and thermal time constants with T, and 
the forcing with H/T. This yields, for steady flows. 



df] 
dx 



-yv- 



dri 
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u 

'^drag 
V 

'^drag 



du dv\ 

ox ay J Trad 



(16) 
(17) 
(18) 



where all quantities, including Trad and Tdrag, are now nondi- 

mensional. 

In p ioneering investigations, Matsuno' (1966) and ^ Gill 
(|1980|) obtained analytic solutions to Eqs. ([T6])-([T8]) for the 
special case where the drag and radiative time constants are 
equal and drag is neglected from the meridional momentum 
equation ([17]); the latter assumption is called the "longwave 
approximation" because it is valid in the limit where zonal 
length scales greatly exceed meridional ones. On tidally 
locked exoplanets, however, the drag and radiative time scales 
can differ greatly, and the longwave approximation may not 
apply, because the flow exhibits comparable zonal and merid- 
ional scales. We therefore retain the full form of Eqs. ([16])- 

Equations ( 16)-([18]) can be combined to yield a single dif- 



ferential equation forthe meridional velocity v (e.g., Wu et al.[ 
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|2Q0T] ): 
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'^drag 



5v_ 
dx 



V = 




(19) 



If on e se eks separable solu tions, then, as described by |Gill| 
( |198Q| ) and |Wu et aL] ( |2QQl| ), the meridional structure oftEe 
solutions to this equation with finite Tdrag are the parabolic 
cylinder functions iljJy). which are simply Gaussians times 
Hermite polynomials^ 



ipriiy) = exp 



_y 

'2V^ J 



Hi) 



(20) 



where V = (rrad/Tdrag)^^"^ is the fourth root of a Prandtl num- 
ber. 

Our goal is to specify the thermal forcing, S(x^y), and solve 
for the unknowns u, v, and r]. In general, any desired pattern 
of thermal forcing can be represented as 



(21) 



«=0 



For tidally locked exoplanets, we expect this pattern to consist 
of a day-night variation in heating/cooling whose amplitude 
peaks at low latitudes and diminishes near the poles. We take 
the forcing to be symmetric about the equator (appropriate 
for a planet with zero obliquity) and, to keep the mathematics 
tractable, retain solely the term 5*0, corresponding to pattern 
of heating and cooling that is a Gaussian, centered about the 
equator, with a latitudinal half- width of the equatorial Rossby 
radius of deformation modified by frictional and radiative ef- 
fects. While the full solution would require consideration of 
Sn for all ^ > 0, the first term, 5*0, will be the dominant term 
for cases where the deformation radius is similar to a plane- 
tary radius, as is the case on typical hot Jupiters. Consider- 
ation of this term alone will therefore suffice to illustrate the 
qualitative features relevant for inducing an equatorially su- 
perrotating jet on tidally locked exoplanets. 

Appendix [B] describes the solution method of Eqs. (p^- 
( [TS] ) and presents the solution for the specific case whereme 
forcing consists solely of the 5*0 term varying sinusoidally in 
longitude, i.e., S(x,y) = Soe^^^i/joiy), where is a constant. 
Figure |2] shows an example for parameter values typical of a 
hot Jupiter or hot super Earth (zonal wavelengths associated 
with the day-night heating contrast of a planetary circumfer- 
ence and radiative time constants of order 10^ sec). For this 
example, the drag time constant is taken equal to the radia- 
tive time constant. Figure [2^ shows the radiative-equilibrium 
height field and Fig. |2|) presents the steady state height and 
velocity fields. 

The solutions exhibit several important features. Although 
the radiative-equilibrium height field is symmetric in longi- 
tude about the substellar point (Fig. [2^), the actual height 
field deviates significantly from radiative equilibrium and ex- 
hibits considerable dynamical structure (Fig. Two fun- 
damental types of behavior are present. First, at mid-to-high 
latitudes (\y\ ^ 1-3 in the figure), the flow exhibits vortical 



The first few Hermite polynomials are //o(0 '- 
4^2-2, and HsiO = - 12^- 



l,Hi(0 = 2^,H2(0-- 




-2 2 

Eastward distance 




-0.01 0.00 0.01 

Accelerations 

Fig. 2. — An example linear, analytic solution for parameters relevant to 
hot, tidally locked exoplanets. (a) Spatial structure of radiative-equilibrium 
height field, heq (orangescale and contours), (b) Height field (orangescale) 
and horizontal wind velocities (arrows) for the linear, analytic solution forced 
by relaxation to the heq profile shown in panel (a) and with nondimensional 
zonal wavenumber k = 0.5 and nondimensional radiative and drag times 
'Trad = '''drag = 5. For a hot Jupitcr or hot super Earth, these correspond to di- 
mensional zonal wavelengths of a planetary circumference and dimensional 
radiative and drag time constants of ~10^ sec (see Appendix |a}. In (a) and 
(b), the horizontal and vertical axes are dimensionless eastward and north- 
ward distance, respectively; one unit of distance corresponds to a dimensional 
distance of one Rossby deformation radius, (^/gH/ 13)^^^. The x marks the 
longitude along the equator where h reaches a maximum and the eddy zonal 
wind changes sign, (c) Zonal (east-west) accelerations of the zonal-mean 
flow implied by the linear solution. The black and dark blue curves give the 
accelerations due to horizontal and vertical eddy transport (terms II and III, 
respectively, in Eq. |22) . The light- green and cyan curves show friction (term 
IV) and the effect ot tne mean-meridional circulation (term I), respectively. 
The red curve shows the sum of all terms. The numerical values adopt a forc- 
ing amplitude Aheq/H = 1. For this value, the nondimensional peak winds 
approach 0.5, corresponding to speeds of ~1 kmsec"^ on a hot Jupiter. 



behavior. The dayside contains an anticyclone in each hemi- 
sphere, manifesting as a pressure high (i.e., local maximum of 
the height) around which winds flow clockwise in the north- 
ern hemisphere and counterclockwise in the southern hemi- 
sphere; the nightside contains a cyclone in each hemisphere, 
manifesting as a pressure low around which winds flow coun- 
terclockwise in the northern hemisphere and clockwise in the 
southern hemisphere. Second, at low latitudes (\y\ < 1), the 
flows are nearly east-west; they diverge from a point east of 
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the substellar longitude (marked with a cross in Fig. ^) and 
converge toward a point east of the antistellar longitude. 

As discussed by |Gill| (p78Q), these features can be inter- 
preted in terms of forceoT^amped, steady equatorial wave 
modes. The mid- to-high latitude feature described above is 
dynamically analogous to that of an ^ = 1 equatorially trapped 
Rossby wave, which exhibits cyclones and anticyclones — 
alternating in longitude — that peak off the equator (see Mat-, 
|iuno 1966, Fig. 4c for an example of the flowfield in this 
mode). The low-latitude feature discussed above is dynam- 
ically analogous to a superposition of the ^ = 1 Rossby wave 
and the equatorial Kelvin wave, which is a fundamental equa- 
torially trapped wave mode with strong zonal winds but very 
weak meridional winds and whose amplit ude is s ymmetric 
about, and peaks at, th e equator (see, e.g., IHoltonl ( |200 4 ) or 
[Andrews et al.| ( |1987| )). Both of these wave modes exhibit 
winds that are primarily east- west at the equator; in the exam- 
ple shown in Fig.|2] the Kelvin component dominates over the 
n=l Rossby component at the equator. 

Equations ([T6])-([T8]) indicate that this is a problem gov- 
erned by two parameters: the radiative time constant and the 
drag time constant. We now examine how the behavior de- 
pends on their values. Figure [3] shows linear solutions, as pre- 
sented in Appendix |B] for dimensionless radiative time con- 
stants of 1, 10, and 100 (top, middle, and bottom rows, re- 
spectively) and drag time constants of 1, 10 and 100 (left, 
middle, and right columns, respectively). For parameters 
appropriate to hot Jupiters (rotation period of 3 Earth days 
and gH ^4x 10^ m^ sec ^), these dimensionless values cor- 
respond to dimensional time constants of ~3 x 10^, 3 x 10^, 
and 3x10^ sec, respectively (see Appendix [a|). When the ra- 
diative and drag time constants are short (upper left corner of 
Fig. [3]), the maximum and minimum thermal (height) pertur- 
bations lie on the equator and are close to the substellar and 
antistellar points; in this limit, the height field is close to ra- 
diative equilibrium (compare top left of Fig. [3] with Fig.[2^p] 
and distinct Rossby- wave gyres do not appear. When tEe 
radiative and drag time constants have intermediate values 
(middle of Fig. |3]), cyclones and anticyclones become visi- 
ble and — as in Fig.j^ — exhibit height extrema that are phase 
shifted westward orthe extrema in radiative equilibrium. Sim- 
ilarly, the height extrema along the equator become phase 
shifted eastward relative to radiative equilibrium; thickness 
variations along the equator become modest relative to those 
in midlatitudes. When the radiative and drag time constants 
are long (lower right comer of Fig. [3]), the height field be- 
comes dominated by the off-equatorial anticyclones and cy- 
clones, with minimal variation of height at the equator. In the 
limits Trad 00 and Tdrag ^ 00, the solution becomes flat at 
the equator and is symmetric in longitude about the x = axis 
(a point demonstrated explicitly in AppendixjC]); Fig.|3]shows 
that this limit is almost reached even for Trad and Tdrag of 100. 

Much of the behavior in Fig. [3] can be understood in terms 
of the zonal propagation of equatorially trapped Rossby and 
Kelvin modes. Kelvin waves exhibit eastward group prop- 
agation while long-wavelength, equatorially trapped Rossby 
waves exhibit westward group propagation. When Trad and 
Tdrag are very short (upper left corner of Fig.|3]), the damping 
is so strong that the waves are unable to propagate zonally. As 
a result, the height is close to the radiative equilibrium height 

Appendix |b] demonstrates formally that, in the limit of either time con- 
stant going to zero, the height field converges to the radiative- equilibrium 
height field. 



field. When the two time constants have intermediate values, 
the propagation produces an eastward phase shift of the height 
field at the equator (the Kelvin component) and a westward 
phase shift of the height field in the off-equatorial cyclones 
and anticyclones (the Rossby component) — exactly as seen 
in Fig. |2]3 and the middle of Fig. |3] As the two time con- 
stants become very long, the westward phase offset of the off- 
equatorial cyclones and anticyclones achieves maximal val- 
ues of 90° . At the equator, however, the height variations go 
to zero; this is explained by the fact that Coriolis forces are 
zero at the equator, so the linearized force balance is between 
pressure-gradient forces and drag. Weak drag requires weak 
pressure-gradient forces and hence a flat layer at the equator. 

Now, the key point of our paper is that these linear solu- 
tions have major implications for the development of equa- 
torial superrotation on tidally locked exoplanets. As can be 
seen in Figs. |2]3 and |3j the wind vectors exhibit an overall 
tilt from northwest-to- southeast in the northern hemisphere 
and southwest- to-northeast in the southern hemisphere. This 
pattern, which resembles a chevron centered at the equator 
and pointing east, is particularly strong when the radiative 
and drag time constants are short, but occurs in all the cases 
shown. This structure implies that, on average, equatorward 
moving air has f as ter-than- average eastward wind speed while 
poleward moving air has slower-than-average eastward wind 
speed, so that u'v' <0 in the northern hemisphere and u'v' >0 
in the southern hemisphere. As shown schematically in Fig.[T] 
this is exactly the type of pattern that causes a flux of eastward 
eddy momentum to the equator and can induce equatorial su- 
perrotation. Since momentum is being removed from the mid- 
latitudes, one would expect westward zonal-mean flow to de- 
velop there. 

The physical mechanism responsible for producing these 
phase tilts are twofold. First, the differential wave propa- 
gation discussed above: this propagation causes an eastward 
phase shift of the height field in the Kelvin waves and a west- 
ward shift of the height field in the Rossby waves relative 
to the radiative-equilibrium height field. Because the Rossby 
wave lies on the poleward flanks of the Kelvin wave, the result 
is a chevron pattern where the height contours tilt northwest- 
southeast in the northern hemisphere and southwest-northeast 
in the southern hemisphere. To the extent that velocity vectors 
approximately parallel the geopotential contours (as they tend 
to do away from the equator when drag is weak or moderate), 
this will generate tilts in the velocities such that u'v' <0 in the 
northern hemisphere and u'v' > in the southern hemisphere. 

The second mechanism for generating the velocity tilts 
needed for equatorial superrotation is simply the three-way 
force balance between Coriolis, drag, and pressure-gradient 
forces. Because drag acts opposite to the velocity, and Cori- 
olis forces are perpendicular to the velocity, this three-way 
force balance requires the velocities to be rotated clockwise 
of -V?7 in the northern hemisphere and counterclockwise of 
-\/r] in the southern hemisphere. Given the expected day- 
night gradients in r], this balance implies that the velocities 
will tend to tilt northwest-southeast in the northern hemi- 
sphere and southwest-northeast in the northern hemisphere. 
We demonstrate this fact explicitly with an analytic solution 
in the limit of Trad ^ in Appendix |Dj even when the height 
field is nearly in radiative equilibrium and hence exhibits no 
overall phase tilts, the velocities themselves develop tilts such 
that u'v' <0 in the northern hemisphere and u'v' > in the 
southern hemisphere (see Fig. [15]). The calculation in the limit 





Fig. 3. — Analytic solutions of linearized shallow-water equations (Eqs. [T6}fl8j , as presented in Appendix [b] for dimensionless zonal wavenumber k = 0.5 
and dimensionless radiative and drag time constants of 1, 10, and 100, corresponding roughly to dimensional time constants of 6 hours, 3 Earth days, and one 
month for parameters appropriate to hot Jupiters (see AppendixjAj. All of these cases are forced by Newtonian relaxation of the height field toward a distribution 
analogous to that in Fig. [2^. 



'Trad ^ is particularly interesting because, in this limit, there 
is no zonal propagation of the Kelvin and Rossby waves: the 
radiative damping is infinitely strong and the zonal phase shift 
of the height field (relative to radiative equilibrium) is zero. 
This is the dominant mechanism for the velocity tilts in the 
top-left panel of Fig.|3] 

To demonstrate explicitly how superrotation would emerge 
from these standing- wave patterns, we analyze the zonal ac- 
celerations associated with these linear solutions. Decompos- 
ing variables into their zonal means (denoted by overbars) and 
deviations therefrom (denoted with primes) and zonally aver- 
aging the zonal-momentum equation (Eq. |9]) leads to (e.g., 
Thuburn & Lagneau|1999l ) 
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where a is the planetary radius and A* =hA/h denotes the 



thickness- weighted zonal average of any quantity A. Eq. ( [22 ) 
is the shallow-water version of the Transformed Eulerian 
Mean (TEM) momentum equation, analogous to that in the 
isentropic-c oordinate form of the primitive equations (see |An- 
|drews et aT 1987, Section 3.9). On the right-hand side, terms 
I, IT, and III represent accelerations due to (i) momentum ad- 
vection by the mean-meridional circulation, (ii) the conver- 
gence of the meridional flux of zonal eddy momentum, and 
(iii) correlations between the regions of eddy zonal flow and 
eddy mass source (essentially vertical eddy-momentum trans- 
port). Within this term, the quantity Ru is the zonal component 
of R (equal to -Qu/h when 2 > and when Q < 0). Term 
IV is frictional drag. The final term represents the time rate 
of change of the eddy momentum. In the linear limit, all the 
terms on the right side of Eq. ( 22 ) have vanishingly small am- 
plitude and, in this case, the solutions in Figs. [2j|3] represent 
true steady states. At any finite amplitude, however, terms I- 
IV are nonzero and would cause generation of a zonal-mean 
zonal flow. 

Figure [2]: depicts these terms for the example solution pre- 
sented in Fig. [2|). As expected, horizontal convergence of 
eddy momentum, term II, causes a strong eastward acceler- 
ation at the equator and westward acceleration in the mid- 
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latitudes (black curve). On the other hand, the acceleration 
associated with vertical eddy-momentum transport, term III, 
is strong and westward at the equator (blue), implying down- 
ward transport of eddy momentum at the equator. The remain- 
ing terms — the mean-meridional circulation (term I, cyan) 
and mass- weighted friction (term IV, light green) — are small 
at the equator. The two eddy terms partially cancel at the 
equator, but the acceleration due to horizontal eddy momen- 
tum convergences exceeds that due to vertical eddy momen- 
tum convergences, leading to a net eastward acceleration at 
the equator and westward acceleration in midlatitudes (red 
curve). 

Remarkably, despite the wide range of morphologies that 
occur when Trad and Tdrag are varied (Fig. [3]), all the solutions 
exhibit an equatorward flux of eddy momentum and a net east- 
ward acceleration at the equator. This is shown in Fi g HI 
which presents the two eddy acceleration terms from Eq. ( [22] ) 
for each of the cases shown in Fig. [3] These solutions there- 
fore suggest that superrotation at the equator and westward 
mean flow in the midlatitudes should occur at essentially any 
value of the control parameters. 

The patterns of spatial velocity and mass source/sink illu- 
minate the physical origin of the westward equatorial accel- 
eration caused by the vertical eddy exchange. The solutions 
show that the longitudes of zero zonal wind at the equator 
lie east of the mass-source extrema (Fig.j2]3), a fe ature also 
clearly visible in the steady, linear calculations of 'Matsuno" 
fl966, Fig. 9) and GiU (1980, Fig. 1). Because of this shift, 
equatorial mass sources (sinks) occur predominantly in re- 
gions of westward (eastward) eddy zonal flow. On average, 
therefore, the mass sink regions transport air with eastward 
column-integrated eddy momentum out of the layer. The mass 
source regions transport air with no relative zonal momen- 
tum from the quiescent abyssal layer into the upper layer; this 
process conserves the local, column-integrated relative mo- 
mentum v/z of the upper layer. Thus, when zonally averaged, 
vertical exchange at the equator removes momentum from the 
layer, leading to (u'Q' -\-Ru ) < and contributing a westward 
acceleration (blue curve in Fig. [2]:). 

The above argument, however, does not determine which 
of the two eddy terms (II and III in Eq. [22]) dominates. To 
determine which is larger — and hence whether the net equa- 
torial eddy acceleration is eastward or westward — we write 
the zonally averaged zonal momentum equation in the form 



du — u 

— = v'C'+v(/+C) +Ru. 

ot Tdrag 



(23) 



where ( is the relative vorticity. For the case where the forc- 
ing is symmetric about the equator, the solutions are symmet- 
ric about the equator in u and h but antisymmetric about the 
equator in v and (. As 3. result, the meridional velocity and 
relative vorticity are zero at the equator, so the terms v(f + Q 
and vX^ vanish there. Therefore, 



du 
dt 



'^drag 



Sity = 0. 



(24) 



Essentially, at the equator, Ru is the mismatch between the ac- 
celerations caused by horizontal and vertical eddy-momentum 
fluxes. The analytic solutions, which assume u = 0, show that 
u is predominantly westward in regions where 2 > 0, which 
therefore implies that Ru > 0. From Eq. ( |24| ), the net eddy- 
induced acceleration is therefore eastward. This explains, in 



a general way, the sign of the net eddy accelerations at the 
equator in Fig. [4] Of course, once a zonal-mean flow (u ^ 0) 
develops, the magnitude of Ru changes and the friction term 
becomes important in Eq. ([24]); eventually these terms bal- 
ance and allow a steady state to be achieved. We discuss the 



possible steady states in light of this equation in 

We have so far emphasized the spatial patterns of the cir- 
culation, but it is also interesting to examine the magnitudes 
of the velocities predicted by our linear solutions. When 
the day-night difference in the radiative-equilibrium height 
is comparable to the mean value and the radiative time con- 
stant is a few days or less (as expected for the strongly forced 
conditions on hot Jupiters), the winds shown in Fig. [2]3 and 
Fig. [3] reach nondimensional speeds of order unity. For a hot 
Jupiter, with typical g - 20msec~^ and H - 200 km, this cor- 
responds to dimensional speeds of ~2kmsec~^ To within a 
factor of a few, this is similar to the speeds obtained in fully 
nonlinear three-dimensional atmospheric circulation models 
of hot Jupiters ( Showman & Guillot 2002; Cooper & Sh ow- 
'man"2005», 'Showman et al. 2008, 2009jPobbs-Dixon & Lm 



2008; Dobbs-Dixon et al. 2010; Men ou & Rauscher 20091 
Rauscher & Menou 2010; Thrastarso n & Cho||2010| )T For a 
tidally locked. Earth-like planet in the habitable zone of an 
M-dwarf, with g = lOmsec"^, H = 10km, and an Earth-like 
radiative time constant of ~10 days (corresponding to dimen- 
sionless time constants of 10-100), the solutions then yield 
nondimensional speeds of ~ 0.02-0.1. This corresponds to 
dimensional speeds of up to a few tens of msec~^ similar to 
speeds obtained in models of tidally locked terrestrial planets 
dJoshi et al.||1997l IHeng & Vogt||2010l IMerlis & Schneider 

^mr — — 

3.2. Nonlinear solutions 

Next, we relax the small-amplitude and Cartesian con- 
straints to demonstrate how nonlinearity and full spherical 
geometry affect the solutions, and we show how the wave- 
induced accelerations interact with the mean flow to generate 
an equilibrated state exhibiting equatorial superrotation. To 
do so, we solve the fully nonlinear forms of Eqs. ([9])-(ll) 
in global, spherical geometry, using a radiative-equmbnum 
thickness given by 

/Zeq = + A/Zeq cos A cos (/) (25) 

where H is the mean thickness, A/Zeq is the day-night contrast 
in radiative-equilibrium thickness, and the substellar point is 
at longitude 0° and latitude 0°. The planet is assumed to 
be synchronously rotating, so that the pattern of /Zeq(A,0) re- 
mains fixed in time. For concreteness, we adopt planetary 
parameters appropriate to a hot Jupiter, although we expect 
qualitatively similar solutions to apply to super Earths. For a 
typical gravity of 20msec~^ and scale height of 200 km appro- 
priate to hot Jupiters, we might expect gH = 4x 10^ m^ sec~^, 
and we adopt this value for all our runs. (Note that g and 
H do not need to be specified independently.) We also take 
Q = 3.2x 10~^sec~^ and a = 8.2 x 10^ m, corresponding to 
rotation period and planetary radius of 2.3 Earth days and 
1.15 Jupiter radii, respectively, similar to the values for HD 
189733b. 

We reiterate that the equations represent a two layer sys- 
tem with an active layer overlying a quiescent, infinitely deep 
lower layer. Because of coupling between the layers (specifi- 
cally, mass exchange in the presence of heating/cooling), the 
solutions readily reach a steady state for any value of the drag 
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Fig. 4. — Accelerations of the zonal-mean zonal flow for the same solutions as in Fig. [3] Black curve is acceleration due to horizontal eddy-momentum 
convergence (term II in Eq. |22) , and blue curve is acceleration due to vertical eddy-momentum transport (term III). Red curve is sum of terms I, II, III, and IV. 
Numerical values plotted here assume that *So/Trad (which is essentially equivalent to Aheq/H) equals one. 



time constant, including the Hmit where drag is excluded en- 
tirely in the upper layer (rdrag ^ oo). This in fact is a simple 
representation of the situation in many full 3D GCMs of So- 
lar System atmospheres, including Earth, which often have 
strong frictional drag near the surface, little-to-no friction in 
the upper layers, and yet easily reach a steady configuration 
throughout all the model layers. In our case, we find that, 
when drag is strong, the solutions reach steady states in run- 
times < lOrdrag. In the case where drag is turned off, the time 
to reach steady state is determined by the magnitude of mo- 
mentum and energy exchange between the layers (e.g., by the 
magnitude of the R term), and is generally < 10rrad|^/A/z|, 
where \Ah\/H is a characteristic value of the fractional height 
variations in the active layer. All solutions shown here are 
equilibrated and steady. 

We solve Eqs. ([9|)-(pT]) using the Spectral Transform Shal- 
low Water Model (STSWM) of H ack & Jakob (1992 ). Rather 
than integrating the equations for u and v, the code solves the 
momentum equations in a vorticity-divergence form. The ini- 
tial condition is a flat layer of geopotential gH at rest; the 
equations are integrated using a spectral truncation of T170, 
corresponding to a resolution of 0.7° in longitude and lati- 



tude (i.e., a global grid of 512 x 256 in longitude and lati- 
tude). A hyperviscosity is applied to each of the dynami- 
cal variables to maintain numerical stability. The code adopts 
the leapfrog timestepping scheme and applies an Asselin filter 
at each timestep to suppress the computational mode. These 
methods are standard practice; for further details, the reader 
is referred to |Hack & Jakobl ( [T992l ). 

To facilitate comparison with the analytic theory in §3.1| 
we first describe the solutions at very low amplitude where 
the behavior is linear. Figure [5] shows the geopotential (i.e., 
gh) for equal radiative and drag time constants of 0.1, 1, and 
10 (Earth) days, respectively. Qualitatively, the numerical so- 
lutions in spherical geometry bear a striking resemblance to 
the analytic solutions on a plane. At time constants of a 
fraction of day, the geopotential maxima occur on the equa- 
tor, and for time constants of 0.1 day (a), the geopotential 
resembles the radiative-equilibrium solution, with wind flow- 
ing from the substellar point to the antistellar point. Longer 
time constants (1 day, panel b) allow zonal energy propaga- 
tion of the Kelvin and Rossby waves, leading to an eastward 
phase shift of the geopotential at the equator and a westward 
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Fig. 5. — Geopotential gh (orangescale and contours) for the equilibrated 
(steady-state) solutions to the shallow- water equations (Eqs. [9}flT) in full 
spherical geometry assuming equal radiative and drag time constants of 0.1, 
1, and 10 days in (a), (h), and (c), respectively. White is thick and dark is 
thin. Although the equations solved are fully nonlinear, the forcing ampli- 
tude is small here (Aheq/H = 0.01) so that the solutions in these panels are 
essentially linear. 



phase shift at high latitudes (~40-90°). The result is con- 
tours of geopotential that develop northwest- southeast tilts in 
the northern hemisphere and southwest-northeast tilts in the 
southern hemisphere. When the time constants are long (10 
days, panel c) off-equatorial cyclones and anticyclones domi- 
nate the geopotential, with only weak geopotential variations 
along the equator. These vortices are oval in shape, exhibiting 
no overall phase tilt, though the regions close to the equa- 
tor do develop phase tilts (westward-poleward to easward- 
equatorward). The momentum fluxes cause a prograde eddy 
acceleration (and superrotation) at the equator for all these 
cases. All of these features are also shared by the analytic 
solutions (Fig. [3]). 

We now explore how the solutions change when the ra- 
diative and frictional time scales are different. In the lin- 
ear limit, the latitudinal width of the region exhibiting pro- 
grade phase tilts (i.e., north west-to- southeast in the northern 
hemisphere and southwest-to-northeast in the southern hemi- 
sphere) contracts toward the equator when the drag time con- 



stant greatly exceeds the radiative time constant. This is 
illustrated in Fig. [6] which shows the equilibrated (steady- 
state) solutions for Trad = Iday and Tdrag/Trad = 1 (top), 10 
(middle), and infinite (bottom). When the time constants 
are equal, the entire northern (southern) hemisphere exhibits 
northwest-to-southeast (southwest-to-northeast) phase tilts. 
When Tdrag/^rad = 10, thcsc phasc tilts are confined within 
~20° latitude of the equator, and for Tdrag ^ oo, the width 
shrinks toward z ero. This behavior is explained by the an- 
alytic theory in |3.1| As shown in Eq. ( [2Q| ), the parabolic 
cylinder functions comprising the latitudinaTstructure exhibit 
a characteristic latitudinal width of L(rrad/Tdrag)^^'^. where 
L = (y/gH / py/'^ is the equatorial Rossby deformation ra- 
dius; these functions thus collapse toward the equator as 
Tdrag /^rad bccomcs infinite p] Poleward of this region, the so- 
lutions exhibit phase tilts oime opposite direction (northeast- 
to-southwest in the northern hemisphere and southeast-to- 
northwest in the southern hemisphere). Appendix |C] gives the 
explanation for this reversal in phase tilts; the low- amplitude, 
full spherical numerical solutions at Tdrag /^rad ^ 1 strongly 
resemble analytic solutions in the absence of drag, presented 
in Appendix [C] (compare Figs. [6]: and [14]). 

Nonlinearity alters the solutions in several important ways, 
which we illustrate in Fig. [7] showing a sequence of solu- 
tions for Trad = 0.1 day, Tdrag = 10 day s and, from top to bot- 
tom, amplitudes A/zeq/// of 0.01, 0.1, and 0.5, respectively. 
This choice of time constants is representative of the regime 
of strong radiative forcing and weak drag that may be ap- 
propriate to typical hot Jupiters. Increasing the forcing am- 
plitude (i.e., increasing Ah^q/H while holding Trad and Tdrag 
constant) of course leads to increased wind speeds and day- 
night geopotential variations; for the parameters in Fig. [Tj 
the zonal-mean zonal wind speed at the equator ranges from 
~10msec~^ at the lowest amplitude to almost 1000m sec~^ 
for the highest amplitude shown. Moreover, beyond a critical 
value of AhQq/H (depending on the values of Trad and Tdrag), 
the solutions begin to deviate qualitatively from the linear so- 
lutions. 

First, nonlinearity allows greater geopotential variations to 
occur along the equator, such that at extreme forcing ampli- 
tude the geopotential extrema can in some cases occur along 
the equator when they otherwise would not. In the linear limit, 
the zonal force balance at the equator in the steady state is be- 
tween the pressure-gradient force and drag (cf Eq.^Sj; there- 
fore, when drag is weak, the pressure-gradient force must 
likewise be small, implying that minimal variations of geopo- 
tential occur along the equator. This restriction does not apply 
at higher latitudes (where the Coriolis force can balance the 
pressure-gradient force), so for very weak drag the thickness 
extrema generally occur off the equator (as can be seen in 
the lower right portion of Fig. [3j Fig. |5|? and c; Fig. [6j and 
Fig.[7]2). At large forcing amplitude, however, the momentum 
advection term v • V/z and the R terms become important and 
can balance the pressure-gradient force, allowing significant 
zonal pressure gradients — and hence significant variations in 
thickness — to occur along the equator. For the parameters in 

The numerical solutions show that the region of prograde phase tilts 
does not become precisely zero as Tdrag becomes infinite because of the R 
term. As shown in Eq. ( [T2) , R plays a role analogous to drag, and the effective 
drag time constant (one over the quantity in square brackets in Eq. [T2) has a 
characteristic magnitude r^adh/iheq-h). This suggests that, in the aosence of 
drag, the region of prograde phase tilts exhibits a latitudinal width of order 
L[(heq- h)/h]^^^ . This goes to zero in the limit of zero amplitude but is 
nonzero at any finite amplitude. 
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Fig. 6. — Geopotential gh (orangescale and contours) for the equilibrated 
(steady-state) solutions to the shallow- water equations (Eqs. [9}flT) in full 
spherical geometry illustrating the effect of varying rdrag in thelinear limit. 
All of the cases depicted have r^ad = 1 day. The drag time constant is 1 day, 10 
days, and infinite in the top, middle, and bottom panels, respectively. White 
is thick and dark is thin. As in Fig.p] Ahoq/H = 0.01 here. 

Fig. [7] the thickness variations peak at the equator when the 
forcing ampHtude is sufficiently large (bottom panel). 

Second, at high amplitude, the phase tilts of wind 
and geopotential tend to be from north west-to- southeast 
(southwest-to-northeast) throughout much of the northern 
(southern) hemisphere — as in Fig. ^ and [7]: — even when 
the phase tilts are in the opposite direction at low amplitude 
(as in Fig. [T]^). This effect can be directly attributed to the 
term R in tne momentum equations. As shown in Eq. ([12]), 
R plays a role analogous to drag. When true drag is weaKor 
absent, the effective drag time constant (one over the quantity 
in square brackets in Eq. [T2]) has a characteristic magnitude 
^^radh/ih^q-h). At large iwcing amplitude, h/QiQc^-h) ~ 1, 
and in that case the effective drag time constant is compara- 
ble to Trad- The linear solutions show that prograde phase tilts 
dominate over much of the globe when the radiative and drag 
time constants are comparable, but when the drag time con- 
stant greatly exceeds the radiative time constant, the phase 
tilts are in the opposite direction (see Fig. |6]). In Fig. [Tj 
'^drag/'^rad = 100, but the ratio of the effective avdig time con- 
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Fig. 7. — Geopotential gh (orangescale and contours) for the equilibrated 
solutions to the shallow- water equations (Eqs.[9j[TT} in full spherical geom- 
etry illustrating the effect of nonlinearity. AIT or the cases depicted have 
'H-ad = 0.1 day and rdrag = 10 days. The forcing amplitude is Aheq/H = 0.01 
0.1, and 0.5 in (a), (b), and (c), respectively. White is thick and dark is thin. 
Range of gh values is 3.97-4.03 x lO^m^sec'^, 3.6-4.3 x lO^m^sec'^, and 
2. 1-5.5 X 10^ m^ sec~^ from top to bottom, respectively. The solution in (a) is 
steady; those in (b) and (c) exhibit small-scale eddy generation and time vari- 
ability (most evident in panel b), although the global properties (total kinetic 
energy, potential energy, and equatorial jet speed) are essentially constant in 
time. 

stant to the radiative time constant decreases from top to bot- 
tom and reaches ~1 in the bottom panel, explaining the tran- 
sition in the phase tilts from Fig.|7]2 through [7]:. Through the 
momentum fluxes that accompany these phase tilts, the equa- 
torial jet becomes broader and more dominant with increasing 
nonlinearity. 

The momentum balance in the equatorial jet can achieve 
steady state in two ways. In steady state, Eq. ([24]) becomes 



= Ru 



'^drag 



(at the equator). 



(26) 



As described previously, when the zonal-mean zonal winds 
are weak, the zonal wind u is predominantly westward in re- 
gions where 2 > 0, so that Ru >0 (^3.1 and Fig. [J]). This im- 
plies an eastward eddy acceleration of the zonal-mean zonal 
winds at the equator, which induces equatorial superrotation. 
In the first type of steady state, corresponding to a regime 
of strong friction (short drag time constant), this superrota- 
tion implies a strong westward acceleration due to friction 
(-t/ /rdrag)- Steady state occurs when the zonal-mean equato- 
rial jet becomes strong enough for the friction to balance the 
eastward eddy-induced acceleration at the equator. We call 



Equatorial superrotation on hot exoplanets 



13 



this the "high Prandtl number" regime. In the second type of 
steady state, which we call the "low Prandtl number" regime, 
the friction is sufficiently weak that the term -w/rdrag is unim- 
portant in the momentum balance. Because of the eastward 
eddy acceleration, the zonal-mean zonal winds can build to 
high speed. Once they do, they change the nature of Ru — the 
larger u becomes, the smaller the extent to which w < in the 
region 2 > 0, as necessary for Ru > 0. Eventually, for suf- 
ficiently large u, the quantity Ru goes to zero at the equator. 
The equatorial jet thus achieves steady state. 

Figure [8] shows examples of each of these regimes illus- 
trating how the momentum balance occurs. The left column 
presents an example with strong drag (Trad = ^drag = 1 day) 
and the right column presents an example with weak drag 
(Trad = 0.1 day and Tdrag = oo); both are equilibrated and steady. 
These are high-amplitude cases, so the thickness (top row, a 
and e) exhibits large fractional variations, and the phase tilts 
exhibit an overall trend of north west-to- southeast (southwest- 
to-northeast) in the northern (southern) hemisphere, as ex- 
plained in previous discussion (cf Fig.[^. These tilts indicate 
transport of eddy momentum from mialatitudes to the equa- 
tor; as a result, the zonal-mean zonal winds are eastward at 
the equator and westward in the midlatitudes (b and/). Inter- 
estingly, however, the relative strengths of the equatorial and 
midlatitude jets differ and reflect the range of possible varia- 
tion. Panels (c) and (g) show the terms in the zonal-mean mo- 
mentum equation, just the spherical equivalent of Eq. ([22]): 
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As expected, horizontal convergence of eddy momentum, 
term II, causes a strong eastward acceleration at the equa- 
tor and westward acceleration in midlatitudes (black curves). 
The vertical eddy-momentum transport, term III (dark blue 
curves), causes a westward acceleration at the equator that 
counteracts the eastward acceleration due to horizontal eddy- 
momentum convergence. In the case of strong drag (Fig. [8]:), 
the cancellation is imperfect, leading to a net eddy-induced 
acceleration that is eastward at th e equator — as predicted by 
the linear, analytic theory in |3 1| (compare to Fig. [2]:). A su- 
perrotating equatorial jet therefore emerges and only reaches 
steady state when the jet becomes sufficiently strong that the 
zonal-mean drag on the jet, -J^Vrdrag, balances the eastward 
acceleration at the equator (Fig. [8]:). On the other hand, when 
drag is absent, the superrotation induced by the eddy fluxes 
becomes quite strong (Fig.jSh. This mean flow alters the eddy 
fluxes, causing them to self-adjust to an equilibrium where 
the accelerations at the equator due to horizontal and verti- 
cal momentum fluxes cancel, leading to no net eddy-induced 
acceleration at the equator in steady state (Fig. [8^). 

We find that equatorial superrotation occurs at all forcing 
amplitudes, even arbitrarily small amplitudes where the so- 
lutions behave linearly. This is illustrated in Fig. l9] which 



shows the equilibrated equatorial zonal-mean zonal wind ver- 
sus forcing amplitude for solutions with a range of Trad and 
Tdrag combinations. We therefore conclude that the mecha- 
nism for generating equatorial superrotation described here 
has no inherent threshold. Nevertheless, other processes — 
not included in the shallow- water model — can in some cases 
overwhelm the desire of the day-night forcing to trigger su- 
perrotation, particularly when the day-night forcin g is weak. 
This occurs for exa mple in the ca ses examined by Suarez & 
Duffy (1992) and Saravanan (1993), where superrotation only 
developed for forcing amplitudes exceeding a threshold value. 
In their case, the tropical wave forcing only triggers superrota- 
tion when it attains sufficiently great amplitudes to overcome 
the westward torques provided by midlatitude eddies propa- 
gating into the tropics. These issues are discussed further in 

In all the cases shown in Fig. ^ where Tdrag is finite, the 
zonal-mean speed of equatorial superrotation scales with the 
square of the forcing amplitude when the forcing amplitude is 
sufficiently small. This in fact is the expected low-amplitude 
behavior in the high-Prandtl-number regime described above: 
at low amplitude, the solutions become linear, such that the 
velocities, height perturbations, and mass source/sink scale 
with the forcing amplitude. Because Ru scales as the product 
of the mass source/sink and the velocities at low amplitude, 
it is therefore quadratic in the forcing amplitude. In the fric- 
tional regime, Eq. ([26]) implies that u at the equator is simply 
Trad^M, and therefore u itself is quadratic in the forcing ampli- 
tude. This behavior breaks down when the solutions become 
sufficiently high amplitude, as can be seen in Fig.[9j The low- 
Prandtl-number regime is more complex and can lead to a va- 
riety of scaling behaviors depending on the parameters. 

The flow in the shallow-water models differ from that in 3D 
models in one major respect. In many three-dimensional mod- 
els of hot Jupiters, eastward equatorial flow occurs not only in 
the zonal mean but at all longitudes, at least over some range 
of pressures. In contrast, although the shallow- water mod- 
els described here all exhibit eastward zonal-mean flow at the 
equator, the zonal wind at the equator is always westward over 
some range of longitudes. This can be seen as follows: Ru is 
essentially the mismatch in equatorial zonal acceleration be- 
tween horizontal and vertical eddy-momentum transport and 
in steady state, when w > 0, will be greater than or equal to 
zero. From the definition of Ru, this implies westward equa- 
torial flow at some longitudes. This trait probably arises be- 
cause the meteorologically active atmosphere has here been 
resolved with only one layer overlying a deep interior; in fu- 
ture work, it would be interesting to explore models that rep- 
resent the flow with two or more layers overlying a quiescent 
interior to see whether they can develop equatorial flow that 
is eastward at all longitudes. 

4. THREE-DIMENSIONAL MODEL OF EQUATORIAL 
SUPERROTATION 

Here, we show how the basic mechanism for generat- 
ing equatorial superrotation identified in ^ occurs also in 
three dimensions under realistic conditions. To do so, we 
analyze t he three-dimensional m od el of HD 189733b pre- 
sented in [Showman eFaL] ( [2009] ). [Showman eTaL] ([2009]) 
coupled the dynamical core of the IVEITgcm ( [Adcroft et al. 
2004| ), which solves the primitive equations of meteorology 
in global, spherical geometry, using pressure as a vertical 
coordinate, to the state-of-the-art, non-gray radiative trans- 
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— Fully nonlinear shallow-water solutions with Trad = 1 day, Tdrag = 1 day, and A/zeq/// = 1.25 {left column) and Trad = 0.1 day, no drag (rdrag — ^ cxd), 
;eq/// = 0.2 (ng/zf column), {a, e) Geopotential gh (orangescale and contours, with thick and thin regions in white and black, respectively). Scale runs 
from 1.1-5.4 x 10^ sec~^ in top panel and 3.3-4.6 x 10^ sec~^ in bottom panel. (Z?, f) Zonal-mean zonal wind u, (c, g) Accelerations of the zonal-mean 
zonal wind in the equilibrated state. From Eq. ( |27) , curves are terms I (cyan), II (black), III (dark blue), IV (light green), (d, h) Mass source gQ at the equator 
(orangescale, white being positive and dark being negative values) and winds (arrows). Mass source scale runs from -41.3 to 45.1 m^ sec~^ (left) and -38.1 to 
33.2m2sec-3 (right). 



fer scheme of [Marley & McKay ( 1999| ), which solves the 
multi-stream radiative transfer equations using the correlated- 
k method to treat the wavelength dependence of the opaci- 
ties. This coupled model, dubbed the Substellar and Plane- 
tary Atmospheric Circulation and Radiation (SPARC) model, 
is to date the only GCM to include realistic radiative transfer 
for hot Jupiters. Gaseous opacities were calculated assum- 
ing local chemical equilibrium for a specified atmospheric 
metallicity, assuming rainout of any condens ates (i.e. ignor- 
ing cloud opacity). [Showman et al.| ( |2009| ) presented syn- 



chronously rotating models of HD 189733b with one, five, 
and ten times solar metallicity and of HD 209458b with solar 
metallicity, along with several models with non- synchronous 
rotation. Their HD 189733b models in particular c ompare fa- 
vorably with a variety of ob servational constraints ( [Showman 



et al.|20Q9 Agol et al. 2010)), and here we focus on their solar 
metallicity, synchronously rotating HD 189733b case. This 
model adopts planetary radius and gravity of 8.24 x 10^ m and 
9.36msec~^. The rotation rate is 3.3 x 10~^sec, correspond- 
ing to a rotation period of 2.2 Earth days. 
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Fig. 9. — Equilibrated (steady-state) equatorial jet speed from the nonlinear 
shallow- water solutions for a variety of forcing amplitudes, radiative time 
constants, and drag time constants. The jet is eastward (i.e., superrotating) in 
all cases. In most cases, the zonal-mean jet speed scales as forcing amplitude 
squared at low amplitude but exhibits a flatter dependence at high amplitude. 
Here, forcing amplitude is defined as gA/zeq/rrad- 

Figure [TO] shows the velocity and temperature structure at 
the 30-mDar level during the spin-up phase of this model — 
after the forcing has had sufficient time to trigger a global 
wave response but before the equatorial jet has spun up to 
high speed. The velocity pattern in the three-dimensional 
model (Fig. [TO]) strongly resembles the standing Kelvin and 
Rossby-wave pattern described in ^ The flow clearly ex- 
hibits the east-west divergence along the equator, emanating 
from a point near the substellar longitude, identified in ^ as 
the standing Kelvin- wave response. The longitude of peak di- 
vergence (i.e., the longitude at the equator where the zonal 
velocity changes sign) lies east of the substellar longitude, 
as expected from the analytic theory and nonlinear shallow- 
water runs in ^ Moreover, the flow exhibits the broad gyres 
in each hemispnere, anticyclonic on the dayside and cyclonic 
on the nightside, identified in ^ as the standing Rossby- 
wave response. As predicted analytically, the velocities in 
these gyres exhibit a northwest-to-southeast (southwest-to- 
northeast) phase tilt in the northern (southern) hemisphere. 
These phase tilts imply that u'v' is negative in the northern 
hemisphere and positive in the southern hemisphere. Eddy 
momentum therefore fluxes from the midlatitudes to the equa- 
tor, and it is this flux that produces the superrotating equato- 
rial jet (see Fig.[T]). The overall qualitative resemblance to the 
analytic calculation in Fig.[2]is striking. 

As in the shallow- water solutions, the three-dimensional 
models exhibit a net downward eddy momentum flux at the 
equator throughout the upper atmosphere where the radia- 
tive heating/cooling is strong. This momentum flux results 
from the fact that, at the equator, (i) the Matsuno- Gill- type 
standing-wave patterns lead to net zonal eddy velocities that 
are predominantly westward on the dayside and eastward on 
the nightside (see Fig. fTQ|, and (ii) net radiative heating oc- 
curs on much of the dayside, leading to net upward veloc- 
ities, whereas net radiative cooling occurs on the nightside, 
leading to net downward velocities. Thus, at the equator, up- 
ward velocities tend to be correlated with westward eddy ve- 
locities and vice versa. This transports eastward momentum 



downward and causes a westward acceleration at the equa- 
tor throughout the upper atmosphere, which counteracts the 
eastward equatorial acceleration caused by latitudinal eddy- 
momentum transport — just as predicted by the analytic and 
numerical shallow- water solutions in ^(see Figs.|2[|4j and 

To quantify the accelerations resulting from these momen- 
tum fluxes, we consider the Eulerian-mean zonal-momentum 
equation in pressure coordinates. By expanding the dynam- 
ical variables into zonal-mean and deviation (eddy) compo- 
nents, and zonally averaging the zonal-momentum equation 
and adopting pressure as the vertical coordinate, we obtaiif^ 
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On the righthand side, the terms describe the meridional 
momentum advection by the zonal-mean circulation, verti- 
cal momentum advection by the zonal-mean circulation, the 
meridional eddy-momentum convergence, the vertical eddy- 
momentum convergence, and friction (represented generically 
by X), respectively. At the equator, the Coriolis term is zero. 
Because of the approximate symmetry of the flow about the 
equator, v and the meridional gradient of u are small there, 
so the mean-meridional advection term is small at the equa- 
tor. The mean vertical-advection term also tends to be weak 
for the flow considered here, and the net zonal acceleration 
at the equator is then determined primarily by a competition 
between the horizontal and vertical eddy-momentum conver- 
gence terms (the analogs of terms II and III in Eq. ( [27] ) for the 
shallow- water system). 

Figure [TT] shows the time evolution of the zonal-mean zonal 
wind aname two eddy acceleration terms at the equator for 
the solar-me tallicity model of HD 189733b from Showman 
et al. ( |2Q09| ). These are vertical averages through the top por- 
tion of the atmosphere where the radiative heating/cooling is 
strong. The zonal-mean zonal wind accelerates rapidly from 
the initial rest state and approaches an equilibrium within 
~100 days (top). As expected, the acceleration due to hor- 
izontal eddy transport is eastward, while that due to vertical 
eddy transport is westward (bottom). Moreover, as suggested 
by the linear and nonlinear shallow- water calculations, the 
magnitude of the horizontal momentum convergence exceeds 
that of the vertical momentum convergence during spin-up, so 
the net acceleration is eastward at early times. A superrotat- 
ing equatorial jet therefore develops. As the jet speed builds, 
the two acceleration terms weaken significantly, and the ra- 
tio of their magnitudes approaches one. As a result, the net 
acceleration drops to zero, allowing the jet to equilibrate to 
a constant speed (top). This model is in the same regime as 
the shallow-water calculation presented in the right column of 
Fig-IU 

The weakening in time of the eddy accelerations seen in 
Fig. [TT] indicates that the mean flow, once it forms, exerts a 
back-reaction on the eddies that alters their structure. The na- 



ture of these changes are illustrated in Fig. 12 The top panel 
shows the temperatures and winds at 30 rnFar pressure after 
the flow at this altitude has become steady; the superrotating 

An equation analogous to this, except using l og-pressure rather than 
pressure itself as the coordinate, can be found in [Andrews et al.| jl987| 
Eq. 3.3.2a). 
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Fi g. 10. — Temperature (c olorscale, in K) and winds (arrows) during the spin-up phase of the three-dimensional, solar-metallicity model of HD 189733b 
from [Showman et al.'f 2009V This shows the state at the 30-mbar pressure level — near the mean infrared photosphere — at a time of one Earth day, after the 
day-night forcing has generated a global wave response but before strong zonal-mean jets have developed. Substellar point is at 0° latitude, 0° longitude. Note 
the development of a "Matsuno-GiH"-type pattern, which leads to an equatorward flux of eddy momentum that pumps the equatorial jet. 
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Fig. 11. — Time evolution of the upper- atmospheric ci rculation in the three- 
dimensional model of HD 189733b by Showman et al.'(^09 ). {Top:) Zonal- 
mean zonal wind at the equator. {Bottom:) Zonal accelerations due to lati- 
tudinal eddy-momentum convergence, -(acos^ (f))~^d(u'v' cos^ 4>)/d(l) {black 
curve) and vertical eddy-momentum convergence, -d{u'uj')/dp {blue curve). 
All quantities are shown in the upper atmosphere, averaged vertically from 
30 mbar to the top of the model. 

equatorial jet, eastward offset of the hottest region from the 
substellar point, and ot her features are evident as detailed in 
[Showman et al. ( 2009 1. The bottom panel depicts the eddy 
temperature diwdeddy winds for the same pressure and time — 
that is, T' in colorscale and iu' ,v') as arrows. Several features 
are similar to those in Fig.fTOj the eddy flow near the equator 
is approximately zonal anoexhibits a Kelvin- wave-like char- 
acter, with predominantly eastward flow at some longitudes 
and westward flow at others; the midlatitudes contain broad 
Rossby-wave gyres in each hemisphere, anticy clonic on the 
day side and cyclonic on the nightside. Interestingly, how- 



ever, the Kelvin-wave structure is shifted eastward, and the 
midlatitude velocity structure differs significantly, relative to 
that with weak mean flow (compare Fig. 12 3 to Figs.j2]3 and 
[T^. From longitudes of about -75° to +4(r , the midlatitude 
velocity structure induces equatorward momentum flux (i.e., 
u'v' negative in the northern hemisphere and positive in the 
southern hemisphere), but at longitudes ~ 80 to 150° the flux 
is reversed (i.e., u'v' positive in the northern hemisphere and 
negative in the southern hemisphere). Due to this cancella- 
tion, the magnitude of the zonally averaged flux u'v' is signif- 
icantly weaker in the equilibrated state than during the spin-up 
phase, when the signs of the midlatitude u'v' add coherently 
at most longitudes (see Fig.[TO]). 

The latitudinal pattern of zonal-mean eddy accelerations in 
the upper atmosphere of the 3D model, shown in Fig. [131 ex- 



hibit a strong relationship to those from the shallow-water 
calculations. The comparison is most apt to shallow- water 
calculations with short radiative time constant (Trad ~ 0.1-1 
day), weak frictional drag (rdrag ^ oo), and large amplitude, 
as depicted for example in the right column of Fig. [8 In the 
3D model (Fig. [13]), the acceleration due to horizontal eddy- 
momentum convergence is eastward at the equator and west- 
ward in midlatitudes, in agreement with analytic theory and 
non-linear shallow-water solutions (compare black curve in 
T3]to Figs. |4]and[8^). Near the poles, the situation is 
Poleward of '^50° latitude, the acceleration 

In the 



Fig 



more complex. 

in the upper atmosphere of the 3D run is eastward, 
shallow- water solutions, the pattern of eddy accelerations at 
high latitude depend on Trad, Tdrag, and the forcing amplitude 
(e.g., compare Figs. [8]: and [8^ poleward of ~60° latitude), 
but for the shallow-water cases most relevant to the 3D run 
shown here — such as the right column in Fig. [8] — the acceler- 
ation due to horizontal eddy convergence becomes eastward 
at high latitudes (Fig. [8^), like that in the 3D run. In the 
3D model, the acceleration due to vertical eddy-momentum 
convergence (blue curve in Fig. [13]) is westward at the equa- 
tor and eastward in the midlatituaes, again like that arising 
in the shallow- water solutions, although a significant differ- 
ence is that the midlatitude eastward acceleration is weak in 
the shallow-water runs but strong in the three-dimensional run 
(relative to the magnitude of acceleration at the equator). 
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Fig. 12. — Top: Temperature (colorscale, in K) and winds (arrows) for the solar-metallicity model of HD 189733b from | Showman et ar]j2009) . Bottom: Eddy 
temperature T' (colorscale, K) and eddy winds (M^ vO (arrows) for the same model. Both are shown at 30 mbar pressure, near the mean mtrared photosphere, 
after the winds at these levels have reached steady state. 
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Fig. 13. — Zonal-mean zonal accelerations from the solar-metallicity 
model of HD 189733b from Showman et al. ( 2009 1. Black and blue show 
accelerations due to latitudinal and vertical convergence of eddy momentum, 
respectively. Both are vertically averaged from 30 mbar to the top of the 
model at 0.2 mbar and temporally averaged from 1 1 to 870 days. 



Of course, the standing eddy patterns and resulting zonal- 
wind accelerations in 3D models depend on the strength of 
radiative heating/cooling and drag, just as they do in the 
shallow- water models. For example, in the shallow- water so- 
lutions, the Kelvin- wave structure and Rossby gyres are spa- 



tially distinct when the radiative and/or drag time constants 
are long and the forcing amplitude is small but not when the 
time constants are short or the forcing amplitude is large (e.g., 
compare the upper left versus lower right of Fig. [3] the top 
versus the bottom of Fig. [5] and the top versus the bottom 
of Fig. |7|). The 3D models shown here lie at an intermedi- 
ate position along this continuum, with Rossby and Kelvin- 
wave structures that are visibly distinct, analogous for exam- 
ple to the shallow- water case in Fig.[2]D. 3D models with very 
strong heating rates, however, seem to exhibit eddy patterns 
lacking distinct Rossby wave gyres, more analogous to the 
top-left case in Fig. |3] and the top case in Fig. |5] Examples 
of models in this regime include t he topmost part of the at- 
mosphere in the models of 'Cooper & Showman (2005 , 2006), 
[&)skinen et al. (2007), Rauscher & Menou (2010) and the 
HD 209458b model of Showman et al. (2009). In contrast, 
cases in the literature with more modest heating rates tend to 
exhibit distinct standing Rossby and Kelvin-wave structu res; 
examples include Showman & Guillot ( 2002, Fig. 5), Heng &| 
Vogt ( 2010, Figs. 1 and 12), the lower portion of some of tEe 
models of K oskinen et al. (2007t Fig. 3b), and several of the 
runs in Thr astarson & Cho (2010), which exhibit a planetary- 
scale cyclone and anticyclone m each hemisphere. 

Despite differences of detail, the overall broad similarities 
described here between the 3D and shallow- water models ar- 
gues strongly that the mechanism for equatorial jet mainte- 
nance that we have identified occurs in both the shallow- water 
and 3D models. 
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5. DISCUSSION 

The development of an eastward equatorial jet — that is, 
equatorial superrotation — is a common feature emerging 
from three-dimensional models of synchronously rotating hot 
Ju piters and extrasolar terrestrial plan ets (Showman & Guil- 
lotll2002l ICooper & Showman""2005" '2006; Showman et aE 
2008, 2009; Dobbs-Dixon & Lin 2008; Menou & Rauscher 
2009; Rauscher & Menou 2010; Pema et al. 2010; Heng et al. 
5010; Jos hi et al. J997; Me rHs & Schneider 2010; Heng 
Vogt 2010|). [Showman & Guillot (2002|^rst pointed out 
that, when the radiative and advective time constants are sim- 
ilar, this superrotation causes an eastward displacement of the 
hottest regions from the substellar point — a phenomenon dis- 
covered on HD 189733b five years later (Knutson et al.|2007j 
|2009 ). Despite its relevance, however, the dynamial mecha- 
nisms responsible for generating the equatorial superrotation 
on tidally locked exoplanets have not been previously identi- 
fied. 

Here, we have shown that the equatorial superrotating jet 
results from an interaction of the mean flow with standing, 
planetary- scale Rossby and Kelvin waves generated by the 
day-night thermal forcing. The strong longitudinal varia- 
tions in radiative heating — namely intense dayside heating 
and nightside cooling — trigger the formation of standing, 
planetary- scale equatorial Rossby and Kelvin waves; this is 
essentially a linear response when wind speeds are modest, 
although nonlinearities affect the wave structure at high am- 
plitude. The Kelvin waves straddle the equator while the 
Rossby waves lie on their poleward flanks. As a result 
of the differential zonal propagation — Kelvin waves prop- 
agating to the east and long-wavelength Rossby waves to 
the west — as well as the multi-way force balance between 
pressure-gradient, Coriolis, advective, and drag forces, the 
velocities develop tilts that resemble an eastward-pointing 
chevron centered at the equator. These velocity tilts pump 
eastward momentum from high latitudes to the equator, 
thereby inducing equatorial superrotation. In steady state, 
the zonal-mean equatorial jet speed near the photosphere 
is determined by a balance between this eastward, wave- 
induced acceleration and westward equatorial acceleration re- 
sulting from vertical eddy-momentum transport and/or drag. 
We demonstrated the mechanism in a hierarchy of dynam- 
ical models — including linear, analytic shallow-water mod- 
els, fully nonlinear shallow-water models, and state-of-the- 
art three-dimensional GCMs. For conditions relevant to hot 
Jupiters, such equatorial superrotation occurs over a wide 
range of radiative heating rates and drag time constants. The 
consistency of the picture emerging from this sequence of 
models with widely varying complexity is encouraging and 
suggests that the mechanism is robust. 

The mechanism identified here has several implications: 

• It implies that the equatorial jet results from a direct, es- 
sentially weakly nonlinear interaction between the ther- 
mally forced waves and the mean flow at the planetary 
scale. Eddy-eddy interactions, including the possibility 
of inverse or forward energy cascades or other turbu- 
lent interactions, may occur but are not essential to the 
basic mecha nism. (This is analogous to the situation 
suggested by |Q' Gorman & Schneider] ( |2007| ) for inter- 
action of baroclinic midlatitude eddies with the mean 
flow on Earth.) 

• The wave-mean-flow interaction produces an equato- 



rial jet whose latitudinal width is comparable to that of 
the Rossby waves, namely the equatorial Rossby defor- 
mation radius modified by radiative and frictional ef- 
fects. For conditions typical of synchronously rotating 
hot Jupiters, this length is comparable to a planetary 
radius, explaining the broad scale of the equatorial jet 
obtained in most hot Jupiter models. 

• When acting in isolation, this mechanism for gener- 
ating superrotation has no inherent forcing-amplitude 
threshold; it operates at any forcing amplitude, unlike 
the behavior repor ted in some earlier studies in the 
terrestrial context dSuarez & Duffy||1992l ISaravanan 
[1993). 

• For parameter combinations appropriate to hot exoplan- 
ets, our models yield flows whose hottest regions along 
the equator lie to the east of the substellar point. This 
property results from the eastward group propagation of 
Kelvin waves. The development of a strong mean flow 
(equatorial superrotation) further enhances the offset by 
its eastward advection of the temperature pattern. To- 
gether, these effects can explain the offse ts observed on 
HD 189733b ( [Knutson et alT^QOTll^QQ^l ). 

Despite the ubiquity of superrotation in our models, the 
flow could shift regimes to one with westward zonal-mean 
equatorial flow if westward equatorial accelerations caused 
by other processes — not considered here — become important. 
For example, if baroclinic instabilities occur in midlatitudes, 
they could cause Rossby wave radiation at midlatitudes, po- 
tentially allowing absorption of Rossby waves near the equa- 
tor. This would contribute a westward wave-induced accelera- 
tion near the equator. If this westward acceleration dominates 
over the eastward equatorial acceleration caused by the day- 
night heating contrast, then the net wave-induced acceleration 
at the equator could be westward and equatorial superrota- 
tion would not occur. This could occur on hot Jupiters if the 
planetary rotation rate is sufficiently fast and heating rates are 
sufficiently low. Note that baroclinic instabilities cannot oc- 
cur in a one-layer model such as the shallow-water model, 
helping to explain why no such transitions to westward zonal- 
mean equatorial flow were seen in the shallow-water models 
presented here. In some cases, Hadley cells may also force 
the circulation into a regime of westward equatorial flow, par- 
ticularly if the planetary obliquity is non-zero; this is the case 
in Earth's troposphere. On the other hand, if the planetary ro- 
tation is sufficiently slow, additional mechanisms for generat- 
ing equatorial superrotati on become possible, even when the 
forcing is ax isymmetric | Del Genio & Zhou||1996 Mitchell 
& Vallis|2Q10| ). Exploring the combinations of orbital semi- 
major axes, stellar fluxes, planetary rotation rates, and atmo- 
spheric compositions under which such regime transitions can 
occur in 3D models is an important goal for future work, as 
such transitions could have major implications for visible and 
infrared light curves. 

It is worth mentioning that the presence of a physical 
surface — and entropy gradients on that surface — promote the 
existenc e of barocli nic instabilities in midlatitudes (see for ex- 
ample Vallis 1 120061 Chapter 6), so tidally locked terrestrial 
planets may be more prone than hot Jupiters to experience 
a regime transition to a flow containing midlatitude eastward 
jets in addition to, or instead of, a flow dominated by equato- 
rial superrotation. 
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In the geophysical-fluid-dynamics (GFD) literature, the 
generation of midlatitude eastward zonal jets is often de- 
scribed theoretically using simple barotropic considerations 
analogous to those summarized in ^ involving the propa- 
gation of waves. While such barotropic considerations seem 
to work reasonably well for midlatitude jets, the barotropic 
framework fundamentally fails to explain the emergence of 
equatorial superrotation in 3D models of synchronously rotat- 
ing exoplanets, where the flow is often steady and symmetric 
about the equator. The theory presented here overcomes this 
obstacle and represents an extension of the barotropic frame- 
work to a thermally forced, stratified, ageostrophic flow with 
finite Rossby deformation radius. 

Specifically, the generation of eastward jets is often in- 
terpreted in terms of the meridional propagation of Rossby 
waves away fr om a source region and th eir dissipation at other 
latitudes (e.g., |HeId|[2000t |Wfis||2006| ). As described in ^ 



this interpretation seems to flow naturally from barotropic the- 
ory, in which free Rossby waves, even at the equator, can 
propagate not only in longitude but also in latitude. In con- 
trast, although our work clearly shows how superrotation can 
emerge on tidally locked planets, meridional wave propaga- 
tion plays no obvious role in our theory. Unlike barotropic 
Rossby waves, the baroclinic Rossby waves in our theory 
are equatorially trapped, confined to an equatorial waveguide: 
they ca n pr opagate in longitude and potentially height but not 
latitude Pj Moreover, under conditions appropriate to typical 
hot Jupiters, these waves exhibit meridional scales typically 
stretching from the equator to the pole. It is not at all clear 
that the paradigm of waves propagating from one latitude to 
another applies in this context. Rather, the velocity tilts that 
allow the meridional momentum fluxes to generate superrota- 
tion appear to be explained by the differential zonal — rather 
than meridional — propagation of equatorially trapped Kelvin 
and Rossby waves. Future work may further clarify the issue. 

APPENDIX 
A. NON-DIMENSIONALIZATION 

The nondimensional solutions to Eqs. (|T6])-([T8]) involve three dimensionless parameters — k, Trad, and Tdrag- Here we relate 
these dimensionless parameters to physicaTproperties for exoplanets. For synchronously locked exoplanets, we expect the forcing 
to correspond to a zonal planetary wavenumber 1, implying a dimensional wavenumber of where a is the planetary radius. 
Thus, our dimensionless wavenumber k has the value aT^i^fgH / that is 



k = 



ma 



1/2 



(Al) 



If we assume for purposes of illustration that the mean layer thickness if is a scale height, then gH is just RT where R is the 
specific gas constant and T is the mean temperature of the atmosphere. For a hydrogen atmosphere where R = 3700 J kg~^ 
this gives 



k = 0J5 



1/4 



1/2 



1/2 



lOOOKy V^days^ 

where P is the rotational period and 7?/ is Jupiter's radius. For a carbon dioxide atmosphere where R = 189 Jkg~^ K~\ 



k=1.2 



lOOOK 



1/4 



3 days 



1/2 



1/2 



(A2) 



(A3) 



where R^ is Earth's radius. Thus, relevant values of k for a hot Jupiter and a hot Earth are similar, in the range of ~0.5-2. Hotter 
atmospheres, longer rotation (=orbital) periods, and smaller planetary radii would promote larger values of k. 

Consider now Trad and Tdrag- Their nondimensional values are their dimensional values times (^/gH/S)^^^. Considering Tnondim 
to be either dimenionless time constant (Trad or Tdrag) and Tdim to be its dimensional counterpart, we have 



'Hiondim "~ '^dim 



1/2 



Again equating gH with RT , for a hot Jupiter with a hydrogen atmosphere we obtain 



"Tnondim — 3.6 



/ Tdim \ / Rj\ 

VlO^secy \ a ) 



while for a hot Earth with a CO2 atmosphere 



'^nondim ■ 



Tdii 



10^ sec 



1/2 



1/2 



1000 K 



1000 K 



1/4 



1/4 



3 days 



3 days 



1/2 



1/2 



(A4) 



(A5) 



(A6) 



Although the theory presented here is for steady, forced conditions, this ^j^-^^ baroclinic equatorial Rossby waves tend to be confined to an eq uatorial 
key distinction holds even for freely propagating, unforced waves: barotropic waveguide. For a discussion of the trapping of equatorial waves, seelH^TtoKl 
Rossby waves can generally propagate in latitude-even at the equator- , ^004 pp. 394-400, 429-432) or Andrews et al. (1987 pp. 200-208). ' ' 
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B. ANALYTIC SOLUTIONS FOR GENERAL trad AND tdrag 

Here we present solutions to the nondimensio nal, linearize d s hallow- water equ ations, ([T6])-([T8]), subject to thermal forcing and 

. For notational brevity, we define a = r^^ and 



Gill 



(1980) and 



Wu et al. 



(2001 



drag. We follow the solution method outlined by 
7 = ^md- Defining 

we convert the coupled equations for u, v, and r] (Eqs.[T6}|T8]) to equivalent equations for q, r, and v: 



dr ^dv ^ . 

'dx~ ~dy " v7JV + ya7r-VQ^^(-^,3;) 



[a dq 
yq-yr+ J - — + 



1 dy y J dy 



a dr 



(Bl) 
(B2) 



(B3) 
(B4) 
(B5) 



Gill| ( |198Q| ) and |Wu et al.lpOOl) neglected the drag term in the meridional momentum equation, which is equivalent to dropping 
the term on the right side of Eq. ( |B5| ). However, we retain the full form of Eqs. ( |B3| )-( [B5] ). These equations are separable, and 
we adopt series solutions 



{q,r,v,S} = ^{qn(x),rn(x),Vn(x),Sn(x)}ilJn(y)- 



Recursion relations for the parabolic cylinder functions are 

dy ~ V 

dy ~ V V^' 



(B6) 

(B7) 
(B8) 



Inserting expressions ( |B6| ) into Eqs. ( |B3| )-( [B5] ), using the recursion relations ( |B7| )-( [B8) ), and invoking the orthogonality of the 
parabolic cylinder functions leads to the system 



dqn 
dx 



dVn 

dx 



n > 1 
n > 1 



Equations (B9 ) and ( Bll ) do not apply for ^ = 0, and Eqs. (B3 HB5 ) instead yield for that case 

dqo 
dx 



(B9) 
(BIO) 
(Bll) 

(B12) 
(B13) 



Given a specified longitude and latitude dependence of the forci ng (and hence given 5'„(x) for all n > 0), our goal is to solve for 
qn{x), rn {x), a nd v„(x). To determ ine g p, we simply use Eq. ( |B12| i. For n=\, Eq. ( |B9| > relates q\ and vq to Si. To determine vq, 
use Eq. l |B13| i. Inserting into Eq. ( |B9| i, we obtain 



(B14) 



Obtaining an equation for qn foxn>2 requires full use of Eqs. (B9HB11 ). First, rewrite Eq. (BIO) as an equation for dr^-ildx 
in terms of r„_2, v„_i, and Sn-i. Next, obtain an equation for r„_2 from Eq. ( |B11| ), an d diffe rentiate this expression to obtain an 
equation for drn-i/dx. Inserting these two expressions into the equation derived from (BIO) yields 



dx \aJ dx \aJ 



2{n-l){a^f/\n-i=-V^Sn- 



(B15) 



Equations ( |B9| ) and (|B15[) form two coupled d iffere ntial equations for qn and v„_i in terms of the known coefficients Si. We solve 
from Eq. ( |B9| ) and ir 



for V, 



. insert this into Eq. (B15 ) to obtain a single ordinary differential equation for qn in terms of Sn, dSn/dx, 
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and5'„_2: 



(B16) 



for n>2. With Eqs. ( B12|), (B14 ), and ( B16 ), all possible qn are determined from the specified forcing terms. To determine v„ 
from the qn, we use Eq. ( |B13 ) for ^ = 0, while for ^ > 1, we use Eq. ( |B9| ): 



dq, 



dx 



n+l • 



(B17) 



To determine r„ from the qn and v„, we use Eq. (B 1 1 ): 

rn = 2{n + 2)qn^2 + 2a^'^ Vn^i- (B18) 

valid for all ^ > 0. 

We now specify the forcing and solve for the response. For simplicity, consider a sinusoidal variation of all the variables in 
longitude: 

{qn{x),rn{x),Vn{x),Sn{x)} = {^n, rn.VnJnje'^'' (B19) 

where qn, ^n, and Sn are complex constants and k is the dimensionless zonal wavenumber associated with the day-night 
heating contrast. We take the forcing to be symmetric about the equator (appropriate for a planet with zero obliquity) and, to keep 
the mathematics tractable, assume that the forcing is represented solely by the term So(x), corresponding to pattern of heating 
and cooling that is a Gaussian, centered about the equator, with a half-width (in latitude) of the equatorial Rossby radius of 
deformation modified by frictional and radiative effects: 



S(x,y) = SoilJo(y)e' 



ikx 



(B20) 



While the full solution would require consideration of Sn for all ^ > 0, the first term, 5*0, will be the dominant term for cases 
where the deformation radius is similar to a planetary radius, as is the case on typical hot Jupiters. Consideration of this term 
alone will therefore suffice to illustrate the qualitative features relevant for pumping the equatorial jet on hot Jupiters. 
With these assumptions, Eq. ( B12| ) implies 

^0 - _ , . ,,9 ^0- 



Given that 5*1 = 0, Eq. (B14) implies that qi = 0. Similarly, Eq. ( |B16| ) implies that 



q2- 



So 



All qn, with n>3, equal zero. From Eq. (B13 ), vq = 0, whereas 



vi = 



(ik-\-^/orf) 



(^7)1/4 

which can be expressed in terms of the real and imaginary components of ^2 = ^2reai + ^^2-,^^ , as 



(B21) 



(B22) 



(B23) 



Vl = [(q^7)^^'^^2, 



(^7)1/4 



(a7)i/4 



^2real]- 



All Vn, with n>2, equal zero. From Eq. (B18 ), 



ro = (4 + 2a^/^y^+2ml)^2 



or equivalently 



^0 



= [(4 + 2a^/2v^)42,eai -2a/:^2.^^J + /[(4 + 2a^/2 v^)^2.^^^ +2a/:^2,eai] 



(B24) 



(B25) 



(B26) 



All r„, w ith n> 1 equal zero. This completes the solutions for q, r, and v. The solutions for u and rj are then determined from 
(see Eqs. [BTp2| ) 



u = 



q + r 

q-r 
2v^* 



(B27) 
(B28) 



It is interesting to consider limits of these solutions as 7 and a become infinitely large (implying that the radiative or drag 
time constants go to zero). The solution presented above implies that, in the limit a ^ 00 at constant 7, or in the limit 7 ^ 00 
at constant a, the quantities q^ So/y/y, q2 So/(2a^^^y), and ro So/^/^. Noting that r] = (2y^)~^[^oV^o(}^) + ^2^2(}^) + 
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roi^o(y)W^^, these Hmits imply that r] S^^y'^e^^^, which, by definition, is simply the radiative-equilibrium height field. Thus, 
when either time constant goes to zero, the height field becomes the radiative-equilibrium height field — even if the other time 
constant is non-zero. 



C. ANALYTIC SOLUTIONS IN THE ABSENCE OF DRAG 



The steady, linearized, nondimensional shallow water equations on an equatorial {3 plane (Eqs. 16]- 18) have a particularly 
simple analytic solution in the case when frictional drag is set to zero. In this case, the nondimensional equation governing the 
meridional velocity (Eq. [19]) becomes 

dv dS 



2 

y V 

'Trad 



dx ^ dx' 



(CI) 



This equation differs vastly from Eq. ( [T9] ) because it no longer contains any y derivatives. As before, we seek separable solutions 
that are sinusoids in x. We specify the forcing amplitude S{x^y) = S(y)e^^^ (where S(y) is a specified function that describes the 
latitude dependence of the radiative heating/cooling), and search for solutions u = u(y)e^^^, v = v(y)e^^^, and r] = fj(y)e^^^, u(y), v(y), 
and fjiy) are complex functions of latitude that we seek to determine. Inserting these expressions into Eq. (|C1|), we obtain 



ik 



V^rad y 



V = ikS 



from which we have the solution 



V = S(y)kTrad 



iy^-kr^ady 

/+^Vad' 



The height field is determined by the fact that dr]/dx = yv (Eq. 16 ), which yields 

/ + //^rrad 



fj = S(y)r,^d'- 



We can then determine u using yu = -drj/dy (Eq. 17 ), which yields 



u = S(y)rr^d 



y^ + /j'^^Trad + 2ikr, 



rad 



rad 



dSiy) 
dy 



'^rad 



y^ + iykTr^d 

y ^-ad 



(C2) 



(C3) 



(C4) 



(C5) 



Note that S{y) can be any function; there is no need to decompose the solution into a summation over normal modes (e.g., 
parabolic cylinder functions), as in the case of finite Tdrag. 
The solutions take on a particularly simple form in the limit Trad oo: 



u = 4Si 



y_ 



V = 

2 
~k 



= -S(y) 
i S(y)y^ 
k 

dSiy 
dy k 



(C6) 
(C7) 

(C8) 



It can be seen that, in this limit, the amplitudes of v are in phase with S. The peaks in r] and u are shifted by 90° in longitude 
relative to S(x^y) and are zero at x = 0. The solution in this case is mirror symmetric about the y axis, unlike the case with finite 
time constants. ~ - 22 - nn 

Figure [14J illustrates the solutions for the case S(y) = Se~^ ^ with / = 0.6 and S being a real constant. As in Figs. |2||3j the 
solutions exhibit midlatitude cyclones and anticyclones, with equatorial flow that is zonal ly al igned and diverges from a longitude 
near the substellar point. Because drag is zero, the layer is flat at the equator (see Eq. |C4I which shows that 7/ = for = 0). 
Interestingly, when Trad is finite, the westward phase shift of the Rossby waves is large at low latitudes and approaches zero at high 
latitudes (Fig. [14])). This leads to phase tilts that are southwest-northeast in the northern hemisphere and northwest- southeast in 
the southern hemisphere — opposite to the cases shown in Figs.[2j|3] In the limit Trad ^ 00, these phase tilts disappear and the 
solution strongly resembles that shown in the bottom right corner or Fig. [3] 

Analysis of the continuity equation explains these behaviors. In addition to the mass source/sink caused by the forcing S, two 
processes affect the local layer thickness — radiative damping and horizontal convergence/divergence. Consider how their relative 
strengths depend on latitude. In the absence of drag, the linearized force balance is geostrophic, i.e., yv -dr^ldx andjw = -dr]/dy 
(see Eqs. [T6}|T7] ). This means that, for a given velocity amplitude, the thickness gradients — and hence lateral thickness variations 
themselves — will be larger farther from the equator. As a result, for a given velocity amplitude, the radiative damping (which is 
-7^/rrad) is stronger farther from the equator. On the other hand, in geostrophic balance, the wind flows parallel to contours of 
constant height, which severely limits the horizontal convergence/divergence; convergence can only come about due to variations 
of Coriolis parameter / with latitude, which are stronger near the equator. In geostrophic balance, the horizontal divergence is 
-/3v//, which is just -v/y for the equatorial beta plane considered here. Thus, for a given velocity amplitude, the amplitude of 
horizontal convergence is large near the equator but small at high latitudes. Given these latitude dependences, we thus expect 
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Fig. 14. — Analytic solutions of the linearized shallow-water equations fT6)-fT8) in the limit of zero frictional drag, (a) shows the radiative-equilibrium height 
field, (b) and (c) show the solutions (height field in colorscale and velocity as arrows) for dimensionless T]-ad=10 and oo, respectively. 



that the thickness changes caused by the forcing (S) will predominantly be balanced by radiative relaxation at high latitude but 
horizontal co nvergenc e/divergence at low latitude. 

Equations (|Cl fjC2| ), which represent the zonal derivative of the continuity equation, describe exactly this balance. The three 
terms correspond to local changes in layer thickness due to the forcing S (right side), mass source/sinks due to radiative damping 
(first term on left side), and changes in the layer thickness due to horizontal convergence or divergence of the fluid flow (second 
term on left side). For a given amplitude of v, the radiative damping increases with latitude (scaling with y), while the effect 
of horizontal convergence decreases with latitude (scaling with y~^). The transition occurs aty ^ V^r^ (or y ~ ^J gHkr^^^/ ^ 
using dimensional quantities): radiative relaxation balances the forcing poleward of this latitude, whereas thickness changes due 
to horizontal convergence balance the forcing equatorward of this latitude. 

When radiative relaxation balances the forcing (y > ^/kr^), the height field is in phase (in longitude) with the radiative 
equilibrium height field. When convergence/divergence balances the forcing (y < ^/kr^), the convergence du/dx + dv/dy, and 
therefore v itself, are phase with the radiative-equilibrium height field. Since dr]/dx = yv, this implies that the height field is phase 
shifted westward by 90° relative to the radiative equilibrium height field. In the transition zone (y ~ Vkr^), these arguments 
imply that the phase lines of thickness — and therefore the velocities themselves — tilt from southwest- to-northeast in the northern 
hemisphere and north west- to- southeast in the southern hemisphere. This explains the phase tilts seen in Fig. 143. In the limit 
Trad ^ oo, the wholc domain lies in the region where convergence balances forcing, explaining the lack of pnase tilts of fj in 
Fig. [15. 
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D. ANALYTIC SOLUTION IN LIMIT trad ^ 

There also exist simple analytic solutions to the dimensionless, linearized shallow- water equations ([T6|)-([T8|), for general Tdrag, 
in the limit Trad 0. This limit is particularly relevant for the hottest of tidally locked exoplanets, whicETdue to high temperature 
and/or significant visible- wavelength opacity, are expected to have short radiative time constants. Moreover, as a simplification 
of the full system it provides insights into the dynamical mechanisms. 

In Appendix |B] we showed that, in the limit Trad 0, the height field converges toward the radiative-equilibrium height field. 
Thus, in the momentum equations, we can replace the height with radiative-equilibrium height. The continuity equation involves 
the difference between r] and r^eq, however, and so we must retain r] in that equation. This leads to the system: 

-—^-yV = (Dl) 

Ox '?~drag 

— +yu = (D2) 

Oy Tdrag 

du^dv\r^ (D3) 

dx dyj Trad 



Equations ( |D1| )-(D2) constitute an algebraic system for u and v that can directly be solved to yield 



u = - (D4) 

(D5) 



''"drag 



— +/rdrag 



Given u and v, Eq. ( [D3| ) can then be solved to determine the (tiny) mismatch between r] and r^eq- 

Consider the limits of this solution for extreme values of the drag time constant. As Tdrag 0, Eqs. (Dl )-(D2) imply that 
u = -^^rig^Veq/ 9x and V = -'^drig^V&q/^y- t^is strong-drag limit, the winds simply flow down the pressure gradient. On the 
other hand, in the limit of large Tdrag, then away from th e equ a tor, we obtain u = -y~^dr]Qq/dy and v = y~^dr]Qqdx — implying 
geostrophic balance. Note that the approximate solutions (|D4|)-(|D5|) are singular at the equator in the limit Tdrag ^ oo; in this 



limit, the divergence becomes infinite at the equator, and the assumption that r] = rj^^ used to derive ( [D4| )-( D5 ) no longer holds. 

Figure p3] displays this solution for a drag-dominated case (Tdrag = 0.01, left column) and case where drag is comparable to 
Coriolis forces over much of the domain (Tdrag = 1, right column). The top row shows the winds with the assumed height field, 
and the bottom row shows the winds and the horizontal divergence. As expected, when drag is strong (left column), the air 
flows directly from the substellar point toward the antistellar point and is almost parallel to the pressure gradients. When drag 
is less dominant, however (right column), the solution forms a Matsuno-Gill-type wind pattern which exhibits velocities that tilt 
northwest-southeast in the northern hemisphere and southwest-northeast in the southern hemisphere. As discussed in ^ this 
pattern of velocity tilts would generate equatorial superrotation. 

The mechanism for generating these velocity tilts is, simply, the three-way force balance between pressure-gradient, Coriolis, 
and drag forces. Because drag points in the opposite direction of the velocities, and Coriolis forces are perpendicular to the 
velocities, this three-way balance must give a velocity direction that is rotated clockwise of -Vr] in the northern hemisphere 



and counterclockwise of -V?7 in the southern hemisphere. A visual inspection of Fig. 15 makes clear that, given the expected 
pattern of Vr^eq, this rotation forces the flow pattern to adopt velocity tilts that are northwest- southeast (southwest-northeast) in 
the northern (southern) hemisphere. 
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on the right. 



Charbonneau, D., Knutson, H. A., Barman, T., Allen, L. E., Mayor, M., 

Megeath, S. T., Queloz, D., & Udry, S. 2008, ApJ, 686, 1341 
Cooper, C. S., & Showman, A. P. 2005, ApJ, 629, L45 
— . 2006, ApJ, 649, 1048 

Cowan, N. B., Agol, E., & Charbonneau, D. 2007, MNRAS, 379, 641 
Crossfield, I. J. M., Hansen, B. M. S., Harrington, J., Cho, J., Deming, D., 

Menou, K., & Seager, S. 2010, ApJ, 723, 1436 
Del Genio, A. D., & Zhou, W. 1996, Icarus, 120, 332 
Del Genio, A. D., Zhou, W., & Eichler, T. P. 1993, Icarus, 101, 1 
Dobbs-Dixon, I., Gumming, A., & Lin, D. N. C. 2010, ApJ, 710, 1395 
Dobbs-Dixon, I., & Lin, D. N. C. 2008, ApJ, 673, 513 
Gill, A. E. 1980, Q. J. Roy. Meteor. Soc, 106, 447 

Hack, J. J., & Jakob, R. 1992, Description of a global shallow water model 

based on the spectral transform method. Tech. rep.. National Center for 

Atmospheric Research Technical note NCAR/TN-343+STR, Boulder, CO 
Harrington, J., Hansen, B. M., Luszcz, S. H., Seager, S., Deming, D., 

Menou, K., Cho, J. Y.-K., & Richardson, L. J. 2006, Science, 314, 623 
Held, I. M. 1999, Bernhard Haurwitz Memorial Lecture, American 

Meteorological Society, available at www.gfdl.noaa.gov/isaac 
— . 2000, Paper presented at 2000 Woods Hole Oceanographic Institute 

Geophysical Fluid Dynamics Program, Woods Hole Oceanographic 

Institute, Woods Hole, MA (available at 

http://www. whoi.edu/page.do?pid= 1 3076) 
Held, I. M., & Hou, A. Y. 1980, Journal of Atmospheric Sciences, 37, 515 
Heng, K., Menou, K., & Phillipps, R J. 2010, ArXiv e-prints 
Heng, K., & Vogt, S. S. 2010, ArXiv e-prints 
Hide, R. 1969, J. Atmos. Sci., 26, 841 

Holton, J. R. 2004, An Introduction to Dynamic Meteorology, 4th Ed. 

(Academic Press, San Diego) 
Joshi, M. 2003, Astrobiology, 3, 415 

Joshi, M. M., Haberle, R. M., & Reynolds, R. T. 1997, Icarus, 129, 450 
Knutson, H. A., Charbonneau, D., Allen, L. E., Fortney, J. J., Agol, E., 

Cowan, N. B., Showman, A. R, Cooper, C. S., & Megeath, S. T. 2007, 

Nature, 447, 183 

Knutson, H. A., Charbonneau, D., Cowan, N. B., Fortney, J. J., Showman, 
A. R, Agol, E., Henry, G. W., Everett, M. E., & Allen, L. E. 2009, ApJ, 
690, 822 



Koskinen, T. T., Aylward, A. D., Smith, C. G. A., & Miller, S. 2007, ApJ, 
661,515 

Kraucunas, I., & Hartmann, D. L. 2005, J. Atmos. Sci., 62, 371 
Leger, A., Rouan, D., Schneider, J., Barge, P., Fridlund, M., Samuel, B., 
Ollivier, M., Guenther, E., Deleuil, M., Deeg, H. J., Auvergne, M., 
Alonso, R., Aigrain, S., Alapini, A., Almenara, J. M., Baglin, A., Barbieri, 
M., Bruntt, H., Borde, P., Bouchy, F, Cabrera, J., Catala, C, Carone, L., 
Carpano, S., Csizmadia, S., Dvorak, R., Erikson, A., Ferraz-Mello, S., 
Foing, B., Fressin, F., Gandolfi, D., Gillon, M., Gondoin, P., Grasset, O., 
Guillot, T., Hatzes, A., Hebrard, G., Jorda, L., Lammer, H., Llebaria, A., 
Loeillet, B., Mayor, M., Mazeh, T., Moutou, C, Patzold, M., Pont, F., 
Queloz, D., Rauer, H., Renner, S., Samadi, R., Shporer, A., Sotin, C, 
Tingley, B., Wuchterl, G., Adda, M., Agogu, P., Appourchaux, T., Ballans, 
H., Baron, P., Beaufort, T., Bellenger, R., Berlin, R., Bernardi, P., Blouin, 
D., Baudin, F., Bodin, P., Boisnard, L., Boit, L., Bonneau, F, Borzeix, S., 
Briet, R., Buey, J., Butler, B., Cailleau, D., Cautain, R., Chabaud, R, 
Chaintreuil, S., Chiavassa, F, Costes, V., Cuna Parrho, V., de Oliveira 
Fialho, F., Decaudin, M., Defise, J., Djalal, S., Epstein, G., Exil, G., 
Faure, C, Fenouillet, T., Gaboriaud, A., Gallic, A., Gamet, P., Gavalda, P., 
Grolleau, E., Gruneisen, R., Gueguen, L., Guis, V., Guivarc'h, V., 
Guterman, P., Hallouard, D., Hasiba, J., Heuripeau, R, Huntzinger, G., 
Hustaix, H., Imad, C, Imbert, C, Johlander, B., Jouret, M., Journoud, P., 
Karioty, F, Kerjean, L., Lafaille, V., Lafond, L., Lam-Trong, T., 
Landiech, P., Lapeyrere, V., Larque, T., Laudet, P., Lautier, N., Lecann, 
H., Lefevre, L., Leruyet, B., Levacher, P., Magnan, A., Mazy, E., Mertens, 

F. , Mesnager, J., Meunier, J., Michel, J., Monjoin, W, Naudet, D., 
Nguyen-Kim, K., Orcesi, J., Ottacher, H., Perez, R., Peter, G., Plasson, P., 
Plesseria, J., Pontet, B., Pradines, A., Quentin, C, Reynaud, J., Rolland, 

G. , Rollenhagen, F., Romagnan, R., Russ, N., Schmidt, R., Schwartz, N., 
Sebbag, I., Sedes, G., Smit, H., Steller, M. B., Sunter, W., Surace, C, 
Tello, M., Tiphene, D., Toulouse, P., Ulmer, B., Vandermarcq, O., 
Vergnault, E., Vuillemin, A., & Zanatta, R 2009, A&A, 506, 287 

Lewis, N. K., Showman, A. P., Fortney, J. J., Marley, M. S., Freedman, 

R. S., & Lodders, K. 2010, ApJ, 720, 344 
Li, J., & Goodman, J. 2010, ArXiv e-prints 
Lian, Y, & Showman, A. R 2008, Icarus, 194, 597 
— . 2010, Icarus, 207, 373 



26 



Showman and Polvani 



Lindzen, R. S. 1981, J. Geophys. Res., 86, 9707 
Marley, M. S., & McKay, C. R 1999, Icarus, 138, 268 
Matsuno, T. 1966, J. Meteorol. Soc. Japan, 44, 25 
Menou, K., & Rauscher, E. 2009, ApJ, 700, 887 
Merlis, T. M., & Schneider, T. 2010, ArXiv e-prints 
Mitchell, J. L., & Valhs, G. K. 2010, ArXiv e-prints 
Norton, W. A. 2006, J. Atmos. Sci., 63, 1420 

O'Gorman, R A., & Schneider, T. 2007, Geophys. Res. Lett., 34, 22801 
Pedlosky, J. 1987, Geophysical Fluid Dynamics, 2nd Ed. (Springer- Verlag, 
New York) 

Perna, R., Menou, K., & Rauscher, E. 2010, ApJ, 719, 1421 
Rauscher, E., & Menou, K. 2010, ApJ, 714, 1334 
Saravanan, R. 1993, J. Atmos. Sci., 50, 1211 

Sardeshmukh, P. D., & Hoskins, B. J. 1985, Quarterly Journal of the Royal 

Meteorological Society, 111, 261 
— . 1988, Journal of Atmospheric Sciences, 45, 1228 
Schneider, T., & Liu, J. 2009, J. Atmos. Sci., 66, 579 
Shell, K. M., & Held, L M. 2004, J. Atmos. Sci., 61, 2928 



Showman, A. R, Cooper, C. S., Fortney, J. J., & Marley, M. S. 2008, ApJ, 
682, 559 

Showman, A. P., Fortney, J. J., Lian, Y., Marley, M. S., Freedman, R. S., 

Knutson, H. A., & Charbonneau, D. 2009, ApJ, 699, 564 
Showman, A. R, & Guillot, T. 2002, A&A, 385, 166 
Showman, A. P., & Polvani, L. M. 2010, Geophys. Res. Lett., in press 
Showman, A. R, Y-K. Cho, J., & Menou, K. 2010, ArXiv e-prints 
Suarez, M. J., & Duffy, D. G. 1992, J. Atmos. Sci., 49, 1541 
Swain, M. R., Vasisht, G., Tinetti, G., Bouwman, J., Chen, P., Yung, Y, 

Deming, D., & Deroo, R 2009, ApJ, 690, LI 14 
Thompson, R. O. R. Y. 1971, Journal of Physical Oceanography, 1, 235 
Thrastarson, H. T., & Cho, J. 2010, ApJ, 716, 144 
Thuburn, J., & Lagneau, V. 1999, J. Atmos. Sci., 56, 689 
Vallis, G. K. 2006, Atmospheric and Oceanic Fluid Dynamics: 

Fundamentals and Large-Scale Circulation (Cambridge Univ. Press, 

Cambridge, UK) 
Watkins, C, & Cho, J. 2010, ApJ 

Williams, G. P. 2003a, Journal of Atmospheric Sciences, 60, 2136 

— . 2003b, Journal of Atmospheric Sciences, 60, 1270 

Wu, Z., Sarachik, E. S., & Battisti, D. S. 2001, J. Atmos. Sci., 58, 724 



